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RVAl -TIME PUMP OPTIMIZATION SYSTEM 
Abstract of the Invention 
A system for optimizing pump operation during oil and gas recovery is provided. By 
strategically disposing a plurality of sensors along the production tubing and sucker rod strings, 
pump operation and performance may be monitored real-time. As an important indida of pump 
performance, dynamic fluid level is provided to the operator/end user on a real-time basis. 
Prerequisite to achieving pump optinuzation, dynamic fluid level and other pertinent data are 
analyzed and enables corrections to be made in the pumping system during operation. A 
computer system having suffident inherent and adaptable expertise is provided to interpret pump 
conditions based upon a plurality of variables and parameters to increase or decrease pump 
production and to maintain a dynamic fluid level determined to be optimal or otherwise 
advantageous by the end user. The system is designed with a panoply of configurations to 
accommodate remote administration of many wells by using serial conmiunication and remote 
transmitting devices. 
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RRA1.-TIME PUMP OPTIMIZATION SYSTEM 
Background of the Invention 

This invention relates to down-hole pump operation, and more particularly relates to 
computer systems for optimizing down-hole pump operation during oil and gas recovery. 

It is well known in the art that downhole pumps are commonly used to provide 
supplemental or artificial lifting action to deliver fluids from subsurface formations to the surface 
of a producing well after reservoir pressure has waned to the extent naturally-available energy is 
insufficient for production purposes. Accordingly, it is common practice for downhole pump 
assemblies and associated pump control systems to be used to transport fluids stored in oil and 
gas wells to the earth^s surface. 

For example, in U.S. Patent No. 5,193,983, Escue et al. teach a pump control system 
having a surface monitoring station for sustaining radio communication with downhole motor- 
pump assemblies. A plurality of sensors incorporated into each downhole motor-pump assembly 
send corresponding signals indicative of such variables as temperature and fluid levels for 
monitoring well and pump conditions. The Elscue disclosure elucidates that prior art pump 
control systems typically monitor minimal variables which have been inadequate for effectively 
identifying the panoply of pump malfunctions which have adversely affected well production. 

For effectively operating downhole pumps during oil and gas recovery, it is imperative that 
the dynamic pumping fluid level be known. As is well known to those skilled in the art, dynamic 
pumping fluid level shows the relationship between pumping rate and well productivity, which, in 



turn, indicates to oilfield practitioners actual well performance. Thus, knowledge of dynamic 
pumping fluid level provides insight into well productivity, completion, and reservoir condition. 
Furthermore, comparison of the actual value for a well's dynamic pumping fluid level with 
theoretical pump output capacity provides crucial insight into the condition and performance of 
the well's pump system. 

As will be readily understood those skilled in the art, a pump loses efficiency as its 
various components wear out. During conventional oil well recovery operations, adjusting pump 
output, which is accomplished by manually adjusting the pump's RPM, is prerequisite to 
maintaining maximum well productivity. As will be appreciated by those skilled in the art, with 
no real-time knowledge of well conditions, analysis of well and pump performance is not only 
limited by lack of timely information, but also is time- consuming and time-intensive. 
Unfortunately, without automated real-time analysis, catastrophic reduction of well production is 
the only way to observe changes pertaining to requirements for pump performance. 

Heretofore in the art, a well operator must use a periodic expensive and inconvenient 
sonic fluid level tiesting to obtain dynamic fluid level. It is common practice, however, due to the 
inconvenience and expense of these sonic fluid tests, for well operators to estimate the required 
production rate — and concomitant pump performance — as an attempt to attain maximum well 
productivity. However, there is a serious risk of excessively low fluid levels will seriously damage 
an oil well's productivity and also causing damage to pumping equipment. To avoid such risk 
that potentially disastrous conditions and consequences will occur, operators frequently take 
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conservative actions, thereby correspondingly reducing the probability of achieving maximum oil 
well production. 

As will be appreciated by those skilled in the art, pumps are used in gas v^lls to remove 
fluid from the well bore, thereby relieving back pressure on the formation. Such back pressure 
produced by fluid accumulated downhole, of course, reduces or may even terminate gas well 
production. These pumps are usually placed by operators below the well perforations* and fluid 
level is reduced until gas flow out of the well resumes. Typically, gas flow emanating from a gas 
well will continue until back pressure recurs due to continuing fluid inflow into the v^ll. 

Thus, it should be evident to those skilled in the art that, by automating a downhole 
pump system and providing continuous analysis of well bore and pump conditions, increased gas 
and oil production of a well heretofore unknown in the art may be attained. Additional benefits 
include increased life span of pumps and wells and cost reduction attributable to optimized well 
production and significantly reduced service requirements. Furthermore, the availabilify of such 
an automated computer system would also assure that optimum longevity and performance be 
obtained from pumps used in the progressive cavity pumping applications and the like. Of 
course, the concomitant advantages pertaining to reduced labor costs and increased production 
performance of the oil or gas well are self-evident. 

It is common knowledge among oilfield production cind automation engineers involved in 
progressive cavity pump systems and the like that there is inherent difficulty associated with 
effectively automating such pump systems under a real-time closed-loop optimization protocol. 
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Not only is the prerequisite data collection instrumentation cost-prohibitive, but also the 
information provided by this collected data must be processed and completely analyzed to assure 
that comprehensive optimization is practicable. Such an oilfield automated pump system has 
been heretofore unknown in the art. Indeed, attempting to construct such a system using 
currendy available instrumentation generally necessitates a practitioner setting pressure sensors at 
the pump suction and discharge* This approach for attempting to automate downhole pump 
systems is evidendy expensive and mechanically-elusive because the pump components being 
controlled are located downhole. 

Accordingly, these linnitations and disadvantages of the prior art are overcome with the 
present invention, wherein a computer system is provided that is particularly useful for enabling 
pumping operations during oil and gas oil well recovery to be monitored and analyzed in real- 
time, wherein production may be optimized. The pump qrstem taught by the present invention 
provides oilfield production and automation engineers with a cost-effective comprehensive 
analytical and automation tool. 

Summary of the Invention 
The present invention provides a system for optimizing pump operation during oil and 
gas recovery. As will be hereinafter described in detail, by strategically disposing a plurality of 
sensors along the production tubing and sucker rod strings, and related downhole apparatus, 
pump operation and performance may be monitored real-time. An important indicia of pump 
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performance has been found to be dynamic fluid level. Accordingly, the present invention 
provides the dynamic fluid level to the operator/end user on a real-time basis. 

Prerequisite to adiieving pump optimization, the present invention analyzes dynamic fluid 
level and other pertinent data and affords a convenient means and method for making corrections 
in the field in the pumping system during operation. It will be understood by those skilled in the 
art that the present invention includes a computer system having suflfident inherent and adaptable 
expertise to interpret pump conditions based upon a plurality of variables and parameters to 
increase or decrease pump production to maintain a dynamic fluid level determined to be optimal 
or otherwise advantageous by the end user. The computer system taught by the present invention 
includes an intuitive user data input interface and pump standard operating performance 
databases. The present invention is designed with a panoply of configurations to accommodate 
remote administration of many wells by using serial communication and remote transmitting 
devices. 

According to the teachings of the present invention, a methodology has been discovered 
which enables downhole pump production to be optimized. 

As will be hereinafter described in detail, it is accordingly an object of the present 
invention to provide an integrated mechanical assembly and computer system for achieving real- 
time downhole pump optimization. 

It is also an object of the present invention to provide means and method for monitoring 
and optimizing dynamic operating condition of a dovmhole pumping system. 
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It is also an object of the present invention to provide means and method for monitoring 
and optimizing dynamic fluid level of a downhole pumping system. 

It is another object of the present invention to provide means and method for sustaining a 
knowledge base of raw data indicative of dynamic operating condition of a downhole pumping 
system. 

It is yet another object of the present invention to provide means and method for applying 
a knowledge base of raw data indicative of dynamic o(>erating condition of a downhole pumping 
system to other downhole pumping systems. 

It is a feature and advantage of the present invention that downhole pumping operations 
are monitored and optimized in a manner and with means heretofore unknown in the art. 

These and other objects and features of the present invention will become apparent from 
the following detailed description^ wherein reference is made to illustrative examples and related 
tables and to the figures in the accompanying drawings. 

In the Drawings 

FIG. I depicts a simplified schematic of the preferred embodiment of the present 
invention. 

FIG. 2 depicts a simplified schematic of a portion of the preferred embodiment depicted 
in FIG. 1. 

FIG. 3A depicts a frontal cross-sectional view of the drive head portion of the preferred 
embodiment depicted in FIGS, i and 2. 

-6- 
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FIG. 3B deplete another frontol crcwMCCtional view of the drive head portion of the 
preferred embodiment depicted in FIGS. 1 and 2. 

FIG. 4 depicts a block diagram of the logic and data flow of the computerized expert 
system of the preferred embodiment of the present invention. 

FIG. 5 depict* a block diagram of the PLC General Loop of the present invention. 

FIG- 6 depicts a bloci diapam of the diagnose procedure of the present invention. 

FIG. 7 depicts a simplified plot of load versus time showing pump slippage treatment 

under the present invention. 

FIG. 8 depict* a plot of pump RPM versus measured axial bearing load for a pumping 

system performing under the present invention. 

FIG. 9 depicts a plot of well head pressure versus pump RPM for a pumping system 
perfbrming under the present invention. 

FIG. 1 0 depicts a plot of fluid level versus pump RPM for a pumping system performing 
under the present invention. 

FIG. 1 1 depicts a plot of fluid level versus pump flow rate fbir a pumping system 
performing under the present invention. 

FIG. 12 depicts a plot of pump performance versus feet of head of water. 

A pump optimization system contemplated by the prefenred embodiment of the present 
invention comprises a mechanical assembly including progreswve cavity pump CPCP*^ means» 
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monitoring means to continuously ascertain the real-time pump performance in a well, and an 
expert computer system for analyzing pump performance and for concomitant adjustment of 
pump characteristics so as to optimize well production. As is known to those skilled in the art, a 
PCP means comprises a stator with a steel tube having an elastomer on the interior of the tube, 
and a rotor which turns inside the stator. Lengths or strings of oil field production tubing 
suspend the stator. Lengths or strings of oil field sucker rods suspend the rotor. Such 
production tubing and sucker rods, of course, are tubular elements routinely used in the oil and 
gas exploration and recovery industiy. 

A drive head provides support for the sucker rods, thereby affording a thrust capability, 
and allows for transmission of rotary torque provided by electric, gas engine, or diesel engine 
power. As will be understood by those skilled in the art, the length of each of the corresponding 
production tubing and sucker rod strings varies with the required pump setting depth relative to 
the location of the oil or gas reservoir. Once installed, the production tubing and stator maintain 
a static vertical elevation in the well: the tubing and stator must critically space the sucker rod to 
assure that proper alignment of the rotor in the stator is sustained. As vrill be appreciated by 
those skilled in the art, the stator provides a slop pin below its elastomeric section to provide an 
indication of rotor location during installation. Oil field pulling units or well service "rigs" 
generally conduct installation. These specialized rigs consist of a mast and mechanical draw- 
works roughly similar to the draw-works used by a conventional crane. 

The comprehensive system for real-time optimization of downhole pump pjerformanoe 
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contemplated by the present invention is depicted in the simplified schematics in RGS. 1 and 2. 
More particularly, FIG. 1 shows the surface components comprising the preferred embodiment. 
In a manner known to those skilled in the art, drive head 125 is fixedly attached to brake 120 
and coupling 1 1 5. Coupling 1 1 5 is coupled to gearbox 1 1 0 which, in turn, is attached to drive 
motor 105. As will be appreciated by those skilled in the art, gearbox 1 10 controls the speed of 
motor 1 05 communicated to drive shaft 1 30. Flow line 5 branches firom the production line at 
flow tee 140 wherein surface discharge pressure is monitored by transducer 25. Similarly, 
transducer 20 is disposed on the annular space to monitor casing pressure. Gas may bubble up 
through casing 150 and, of course, in a manner known in the art, such casing may run to a gas 
gathering system or may be vented to the atmosphere away from the well to avoid safety hazards 
and the like. Axial load on drive head 125 is monitored by load cell 180. 

As clearly shown, output from surface discharge pressure transducer 25, casing pressure 
transducer 20, and load cell 180 are electrically communicated to programmable logic controller 
("PLC") 400 or a field personal computer ("PC") containing an integrated PLC, a remote 
control unit C*RCU''). and an advanced control unit ("ACU"), as will be hereinafter described 
in detail. Communication link 85 is interconnected with PLC 400 and variable frequency drive 
80, which controls the speed of motor 105. Power line 75 is in electrical communication with 
drive motor 105 and variable speed frequency drive 80. The present invention also contemplates 
communication between PLC 400 and the like with a remote PC 450 containing an ACU, 
enable by radio link, modem, direct cable connection, etc Depicted therein is production tubing 
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String 145 contained within casing 150 common to the downhole art. 

The real-time pump optimization system of the present invention thus tracks variables 
providing insight into the relationship between pump discharge pressure and downhole thrust at 
the drive head. Establishing the amount of the load carried by the drive head thrust bearing has 
been found to be a key to performing this analysis particularly on a real-time basis. Preferably, a 
resisbve bridge load cell measures this drive head thrust bearing load. Analog outputs from other 
pressure transducers provide surface discharge pressure and casing pressure. Motor current on 
the drive head motor is simultaneously monitored. As will be appreciated by those skilled in the 
art, these analog values are analyzed to establish the dynamic pump operating condition of the 
pumping system, to establish the dynamic fluid level in the well, and to provide raw data for 
analyzing operating condition of the pump system. 

In the preferred embodiment of the present invention, a PLC or industrial-strength PC 
acts as the communication center linking the mechanical drive and the expert computer system. 
As will be understood by those skilled in the art, the PLC is configured to store data and to 
operate the pump system based on basic information downloaded from the expert computer 
system. The PLC also provides the interface for administering pump performance from a remote 
site. It should be evident that, for a non-automated application of die present invention, the 
PLC may be excluded; under such circumstances, of course, the expert system of the present 
invention would preferably interface with data-providing instruments at the drive head direcdy 
from a portable, field-carried computer. As appropriate, in some applications, an industrial PC 
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configured with analog and digital I/O may be used at the field site. It should be apparent that 
the expert system taught by the present invention would be loaded on this PC, wherein the 
complete pump optimization system operates in the field at the local site. Using technology 
known in the art, the PC-contained system or PLC may be administered from a remote location 
via a modem or radio link. 

Referring now to FIGS. 1 and 2, there is depicted a simplified schematic of a preferred 
embodiment of the downhole rotor-stator assembly which comprises a medianical aspect of the 
present invention. As is well known in the art, the term "downhole" is contemplated to mean 
that an encased device is placed below ground within the annular space of an oil or gas reservoir 
or well. In a conventional manner also well known to those skilled in the art, casing extends 
downwards from the wellhead and is perforated at its lower end to enable formation fluid to flow 
therein and then be forced upwards toward the surface. In particular, formahon fluid is 
conducted to the surface by flowing to the surface inside a production tubing string contained 
within the casing. Packing means is generally used to seal the annulus between this casing and 
the production tubing string. Rotor-stator assembly 1 00 of the oilfield pump system taug^it by 
the present invention comprises rotor means 1 60 and stator means 1 70. Rotor means 160 is 
preferably constructed firom a high strength, precision-machined, chrome-plated steel external 
helix. Stator means 1 70 consists of an internal helix whidi is preferably predsion-molded from a 
durable synthetic elastomer. A conventional oil and gas recovery installation incorporates sudi a 
stator means into the production tubing string 145. 
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Plurality of Airerican Petroleum Institute ("API") sucker rods 1 55 are configured to 
suspend rotor means 1 6> wHthin stator means 1 70 and to drive rotor means 1 60 rotationally. It 
will be understood that sucker rods 135 suspend rotor means 160 within conresponding stator 
means 1 70 and drive rofeor means 160 in a rotational direction. That is, as shown in FIG. 2. 
rotor 160 is driven by sucker rod string 153 which is connected at its lower end to rotor 160 and 
extends inside production tiibing 145 up to the surface. Sucker rod string 1 55 is driven in a 
rotary manner by surface drive head 125 that actuates pump means 100. Tubing string 145 
contemplated hereunder secures stator 1 70 as a stationary member of the pump assembly, at a 
fixed subsurface elevation. FIG. 2 depicts the pump assembly at a level in excess of 6,000 feet 
below the surface. As wall be appredated by those skilled in the art, when rotor means 160 and 
stator means 170 are in place, seal cavities are formed. Then, as rotor means 160 turns in 
corresponding stator means I 70» these seal cavities progress in an upwards direction to discharge 
pumped fluid into tubing string 145. 

Progressive ca-vily pump means 165, a basic component of the preferred embodiment, 
consists of single helical wotor I 60 engaging double helical stator 1 70 that is altodied to the 
bottom of tubing string 1 45. Rotor 1 60 is typically attached to sucker rod string 1 55 that is 
suspended and rotatedby surface drive 125. Surface-mounted drives 125 support and rotate 
sucker rod string 155 thereby transferring torque to downhole progressive cavity pump 165. 
Rotary motion is normJIy obtained through the action of a pulley and beU drive system that can 
be either fixed speed <x variable speed. Variable speed, of course, may be either mechanical or 
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electrical. As the rotor turns eccentrically within the stator of a progressive cavity pump 
contemplated by the present invention, a series of sealed cavities fonms and progresses from the 
suction end to the discharge end of the pump in a manner well known in the art. Accordingly, a 
continuous positive displacement flow is engendered having a discharge rale proportional to the 
rotational speed of the rotor and the differential pressure across the progressive cavity pump. 

Located at the surface of the well is drive head 125 that comprises bearings, seals, etc, 
which are required to rotate plurality of API rods 1 55 in a manner known in the art, thus turning 
rotor means 160 within stator means 1 70. Since progressive cavity pump 165 contemplated by 
the present invention is a positive displacement pump, as the rotabonal speed of pump means 
1 65 varies, the pump output varies proportionally. As will be appreciated by practitioners in the 
art, oilfield progressive cavity pump applications range significandy as a function of the setting 
depth of the pump assembly and the pumping rate prerequisite to sustaining the intended fluid 
output. 

Thus, the preferred embodiment of the present invention comprises an assembly of 
mechanical and raw data collection devices. Drive head 125 comprises a support structure which 
holds pump system drive shaft 130, and thrust and radial bearings to provide a conventional 
mechanism to rotate plurality of sucker rods 155 which rotate rotor 160 in stator 1 70. This 
assembly provides isolation of pumped fluid F by the means of packing gland or seal 135 to 
provide manageable disdiarge of fluid F. Drive head 125 contains radial and axial bearings that 
conventionally support load L attributable to plurality of rods 155, fluid column F, and pump 
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1 65. As is well knowa in ihe art, bearings are used to centralize drive head shaft 130. Drive 
head 125 makes it possible to isolate the elements that bear the load in the axial thrust bearings, 
and. of course, the drive head itself This isolation helps measure axial load L by means of 
hydraulic or electronic instrumentation 400 that is placed between the drive head thrust bearing 
and the drive head. is conventional in the art, drive head 125 is configured with conventional 
devices such as coupling means 1 1 5 and gearbox 1 10 to receive various attachment motors 1 05, 
engines and other conunon prime movers. 

Referring now to the frontal cross-sectional view of drive head 125 and assoaated 
components depicted in RGS. 3A and 3B, it will be understood that an important aspect of the 
present invention is the mechanical relationship engendered by the affect of pump discharge 
pressure and downthrust at the drive head relative to the downhole pump system. Indeed, as will 
be hereinafter described, to monitor this mechanical relationship it has been found to be 
advantageous to establisli a value for suction pressure existing downhole. This data is ascertained 
by measuring the amount of the load that is carried by the drive head thrust bearing. As will also 
be understood by those sldlled in the art, die present invention uses a suitable load cell, such as a 
resistive bridge load cell or a hydraulic load cell and an analog pressure transducer to measure 
this load. 

Thus, drive head 125 conventionally includes a bearing chamber for holding bodi radial 
bearings and axial bearings, and a is top cover for sealing the system from environmental 
contaminants and for ontaining the radial bearing. It will be appreciated that the axial bearing 
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bears the weight of fluid, the plurality of rods, and the pump load. The radial bearings centralize 
the drive head shaft and provide radial load capacity to the system. As is well known in the art, 
the top seal means and bottom seal means work in a>mbination to keep grease inside and dirt 
outside. Braking means 120 prevents back spin when power is no longer driving the sucker rod 
rotationally. Stuffing box or shaft packing means 135 provides a medianical fluid seal between 
the atmosphere and fluid filled tubing 145. This assembly thus enables isolation of pumped fluid 
F by means of packing gland or seal 1 35 affording manageable discharge of the fluid. Also 
shown is load cell 180 which is a mechanical load measuring device. 

In the preferred embodiment of the present invention, as shown in the frontal cross- 
sectional views depicted in FIGS. 3A and 3B, load cell 200 comprises a hydraulic or electronic 
load cell that is placed beneadi thrust bearing cone 230 and cup 225. As will be appreciated by 
those skilled in the art, there are two prevalent load cells commonly used in the oilfield art: a 
hydraulic load cell and a strain gauge load cell. Conventional hydraulic load cell 200 depicted 
therein is contained within casting means 245, shown relative to shaft 130 and shaft sleeve 240. 
Load cell 200 comprises piston 220 and a corresponding load cylinder, silicon O-ring seals 2 1 5, 
pressure transducer 205, and a conventional inlet valve. Also shown are purge plug 260, bearing 
clearance spacer 250, oudet means to pressure pump and valve 2 1 0, and load cell receiver means 
255. As will be appreciated by those skilled in the art, hydraulic load cell 200 is disposed below 
the axial bearing and the drive head shaft, and provides support thereto. In a manner well known 
in the art, piston 220 of hydraulic load cell 200 lifts the axial bearing so that load cell piston 
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bears the same load as the ajcial bearing. As should be evident to those skilled in the art, this 
load corresponds to the veight of total hydraulic fluid load plus the weight of the plurality of rods 
submersed in tlie fluid. The pressure transducer registers this force (attributable to the load) on 
the load cell piston. The axial load of the system is known because of the fundamental 
relationship between pressure and force: 

Load, LB =paxPh8 

where LB is the load of syslem due to the sum of weight of fluid column, rods, the pump, etc.; pa 
is the area of the load cell piston; and Phs is the pressure registered by the load cell. 

Now referring specifically to FIG. 3B, in another embodiment of the present invention, load 
cell means 270 comprises an electronic strain gauge load cell comprising button strain gauge 275, 
load cell frame 285, and contact pins 280. According to the present invention, this system is found 
below and in support of the axial bearing and the drive head shaft, wherein strain gauge 275 is 
placed on the edge of load cell frame 285. Two contact pins 280 with the same altitude are placed 
on the same oenterline from shaft 1 30 at 120*" intervals, thereby being disposed ecjuidistantly from 
button strain gauge 275. The strain is positioned 0 or 360** with the first pin at 120* and the 
second pin at 240**. As shodd be evident to those skilled in the art, this placement allows an equal 
load to be distributed amc^rg these three elements. The load cell frame bears the same load as the 
axial bearing. This force Ooad) on the load cell frame is transferred to the button load cell and the 
two contact pins. The button load cell supports one-third of the total axial load. The strain gauge 
load oell registers this load. It will be appreciated by those skilled in the art that an advantage of 
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measuring one-third of the total axial load is the reduction of size and associated cost of the strain 
gauge element imparting physical and economic feasibility to the design. 

Also depicted in FIGS. I -3A is well head discharge pressure transducer 205 which provides 
an analog value to represent the surface discharge pressure. As will be understood by those skilled 
in the art, surface discharge pressure is the pressure required to overcome surface restriction or back- 
pressure existing in the flow line of the gathering system- This surface discharge pressure 
corresponds to a variable dependent upon such factors as well fluid viscosity, number of operating 
wells (disdiarging into the same gathering system), flow line size, flow rate, elevations changes, etc 
Casing pressure transducer 1 20 is shown for providing an analog value representing the pressure on 
the annular space due to gas associated with petroleum or gas production. Variable frequency 
inverter 80 provides the ability to adjust the speed of the pump by changing the finequency of the AC 
voltage supply to the induction motor turning the drive shaft at the drive head, and provides motor 
current feedback to the computer system taught by the present invention. Also shown is drive motor 
1 05 which provides rotary force to drive shaft, rods and pump components. As will be dear to those 
skilled in the art, drive motor 105 can be configured as a direct drive, gear motor chain drive, etc., 
as appropriate for optimum pump performance as contemplated herein. 

Referring now to the block diagram depicted in FIG. 4, the logic and data flow 
characterizing the expert computer system contemplated by the present invention are depicted. As 
will be understood by those skilled in the art, the computer ^tem taught by the present invention 
is comprised of a field instrumentation unit or a field instrumentation marshall ("FIU") 1 0, a 
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programmable control ("PLC) or remote control unit ("RCU**) 400, and an advanced control 
unit C*ACU") 450, Also depicted are variable speed controller 60, PLC historic data unit 415. 
software interface unit 430, simg program/productivity analysis unit 500, and real-time analysis unit 
470. 

More particularly, referring now to FIGS. I -2 and 4, FIU 1 0 receives inputs from load cell 
180, flowiine pressure transducer 25, casing pressure transducer 20, motor temperature transducer 
30, and additional variables 35 as herein described. As also hereinbefore described, surface drive 
head 125 provides the neasurement for the axial load, supported by the axial bearing. FIU 10 
consists of the well surface instruments necessary to transmit all of the well operational data to the 
RCU 400, and is based on industrial instrumentation standards. The RCU 400 performs the basic 
control rules to optimize the progressive cavity pump system as contemplated under the present 
invention to protect the punp operation from extreme conditions and to assure the well's continuous 
operation. The RCU also performs the communication interface between the ACU 450 or 
SCADA System, based oxi standard protocols. The ACU 450 corresponds to the expert system 
contemplated under the present invention, performing the advanced control rules to optimize well 
production, and performs advanced analysis and diagnostics based on the real-time information 
coming from the field. The ACU generates control actions over the well and alarm messages on an 
operator console, followed by a detailed explanation thereof. Representative control actions indude 
adjusting RPMs, adjusting opening percentages for a valve, etc As will be hereinafter described, 
the ACU 450 through itsreal time anaJysis 470 also provides the user interface with the pump sizing 
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and performance modules 500 contemplated under the present invention. 

As will become clear to those skilled in the art, basic rules as contemplated by the present 
invention comprise rules that adjust pump speed to control production rate to obtain an intended 
fluid level and to provide alarms and/or shutdowns to prevent damage to expensive pump system 
components. Motor speed and concomitant pump output 80 is determined via variable control 
device 60, e.g., servomotor control, using control data to/from other system devices 70 as described 
herein including diluent control, local alarms, and other site devices that require automation. On the 
other hand, advanced rules comprise rules generated through the representation of a knowiedge-based 
oilfield experience, wherein production may be maximized, the chanc3e of operational failures may 
be minimized, potential pump failures and well condition may be ascertained, and appropriate 
preventative actions for optimized pump operations may be recommended. 

As will be understood by those skilled in the art, the ACU taught by the present invention 
can use RCU information, whether it is directly integrated thereto or interconnected to a SCADA 
system. Standard interface protocols 430 for stand-alone or network versions such as DDE, 
TCP/IP, etc., are supported. G)nnection options include radio links, dedicated or a non-dedicated 
phone-line modem, or direct connection 90. It will, of course, be understood that the ACU and 
the RCU may be linked in a dedicated industrial PC located in the field. It will also be understood 
that the ACU can network, control, and operate a plurality of pump installations contemplated 
under the present invention. 

Thus, it will be appreciated that the present invention provides a computer system (see FIG. 
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5) for Optimizing downhole pump performance using real-time data provided by inexpensive, 
commonly available transducers and the like that are strategically disposed within the well head 
assembly as hereinbefore described in detail. Generally, those skilled in the art have sou^^it to 
automate real-time closed-loop PCP operation, but have been unable to avoid using expensive 
instruments and the like which have typically been emplaoed downhole. Besides providing 
sufficiently accurate and current data* such an optimization system must also efficiently process and 
analyze such data so that suitable adjustments may be made in the field in real-time or at least near- 
real-time. 

As will be hereinafter described, software aspects of the preferred embodiment of present 
invention have been developed in C++ for Microsoft Windows 3.1, 95, or NT platforms, and is 
compatible with standard TCP/IP, DDE communications protocols, and any application compliant 
with these protocols. It has been designed with a Client/Server arciiitecture. Of course, any other 
suitable implementation language on any other platform in the art is within the teachings of the 
present invention. Artificial Intelligence techniques such as neural networks, fuzzy logic, and genetic 
algorithms are used to perform the inference engine tasb of the present invention. By way of 
example, as will be hereinafter described, a fiizzy logic module manages the task of comparison and 
detection of the variables' conditions. A fuzzy set consisting of four conditions for each variable is 
used to prevent ambiguous conditions from occurring. 

Another aspect of the present invention is the automation analysis and software comprising 
the expert system and related systems. In the prefened embodiment, analytical tool for progressive 
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cavity pump systems design 22 comprises software that develops mathematical models of a 
progressive cavity pump well to decide the well potential. Well potential is contemplated by the 
present invention to connespond to inflow performance ratio (IPR) and the outflow performance ratio 
(OPR). The correct size of the progressive cavity pump equipment to be used for a well, e.g., 
pumps, drive heads, rods, etc, is also ascertained by this aspect of the present invention. This 
analytical tool works on-line, and off-line, and indudes a database that makes saving and retrieving 
actions possible. 

Now focusing on well potential, the nodal analysis module of the preferred embodiment 
allows the generation of both in-flow and out-flow performance ratio curves for the well and well 
completion. It is possible to set and fix the operating point of the well for analysis of the design, 
redesign, or operating parameters. Operating parameters include intended rate or fluid level, 
hydraulic power, mechanical power, electrical power. Studying the well completion is possible using 
both vertical and horizontal multi-phase flow correlations for li^t, medium, heavy, or extra-heavy 
crude oil. By using these analyses, the optimum intended production rate for the well may be 
determined. Additionally, deciding the necessary pressure differential is possible (Delta P) across 
the pump; Delta P defines the hydraulic power required by the pump. As will be hereinafter 
described, determining the suction pressure (fluid level) for optimum well performance is also 
possible under the present invention. Using a panoply of mathematical correlations, the mechanical 
loads managed by the drive head and the motor (or other prime mover) can be predicted. 

In order to accomplish the pumping optimization objectives of the present invention, friction 
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and density calculations are prerequisite for establishing, in turn, accurate calculations for fluid levels 
because of the impact of friction and density upon flow-line or surface-discharge pressure. As will 
be understood by those skilled in the art, friction and density also contribute hydraulic load and 
resultant mechanical load when involved in pumps producing through a tubular pumping system as 
contemplated hereunder. Friction losses are considered mechanical in nature because of the natural 
resultant increase in pressure due to fluid flow through pipe; these values vary with fluid 
characteristics such as viscosity, tubular dimensions, and pump flow rates. Fluid density calculations 
obtain the specific gravity or veight of a column of fluid. Establishing values through calculations 
for the fluid density allows the exp>ert system aspect of the present invention to mathematically 
establish increased hydraulicloads and buoyancy-impact on suspended sucker rods. As will become 
evident to those skilled in the art, these values are required to accurately calculate the impact of 
system variables and to make corrections for accurately enabling fluid level and mechanical 
interpretations to be made. 

Time-related variables used in the computer system taught by the present invention define 
a predictable mechanical relationship. Indeed, these relationships indicate the potential requirement 
for adjustment to the opcra.ting conditions of the pump or well dynamics. Additionally, these time- 
value relationships provide historic data for use to analyze pump and well conditions preferably on 
a real-time basis. As will be appreciated by those skilled in the art» these analyses define event 
characteristics which may then be stored and used by the Expert System to control, optimize, and 
predict required adjustnnents to assurance optimal pump performance. Of course, these events also 
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indicate demands for service and other pump system necessities. Thus, the present invention affords 
the operator the ability to effectively increase well productivity, reduce down-time, and substantially 
improve operating economy. 

The present invention also addresses the equipment sizing issues that are prevalent in the 
oilfield. In particular, pump sizing and elastomer type determinations are made based on the 
pressure and production requirements. The selection is based on the crude oil chemical affinity to 
available elastomers. As examples, elastomers may be selected having reduced gas permeability or 
resistance to aromatic hydrocarbon or increased maximum temperature capacities. Similarly, drive 
head sizing is performed based on the maximum load held by the axial bearing.and mechanical power 
requirements. Mechanical loads are determined by hydraulic requirements of the pump and weight 
of the rods. Motor selection is based on the maximum amount of power required by the pump. Rod 
string sizing is based on depth, fluid viscosity, API gravity, torsional load, etc. 

Another aspect and advantage of the present invention is rules module 23 based on expertise 
for optimized operation. Such a collection of rules preferably comprises a dynamic knowledge base 
which derives from on-field experience to optimize production and diminish failure possibilities. 
Real-time information constituting field variables including axial load, current motor temperature, 
well head pressure, and rpms, torque, etc., is communicated to the software aspect of the present 
invention means of a standard communications platform. As will be hereinafter described, an 
analysis is made on that raw data to decide the appropriate control action to be taken. As will be 
dear to those convereant in the art, such actions include change of pump speed (RPM) or change 
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of Opening percentage of a bj^ass valve, pump shutdown, and system alarms^ etc. It will be 
appreciated that optimization objectives are to achieve continuous, uninterrupted well production; 
to predict and correct malfunction situations; to lower operating costs; to maximize useful life of the 
equipment by protecting them or by modifying operating conditions. 

As will become evident to those skilled in the art, the diagnosis provided by the present 
invention (see FIG. 6) is accomplished by an integrated and iterative process of pattern recognition 
implemented by a panoply of artificial intelligence tools including neural networks, genetic 
algorithms, fuzzy logic, expert systems, etc. Thus, the on-line closed-loop optimization contemplated 
by the present invention depends on certain variables whidi should be automated in a progressive 
cavity pump wiJl such as: system axial load, rotor RPM, motor current, automatic shutdovm and 
alarms, chemiail and diluent injection control, motor winding temperature. As is common in the art, 
production engineers must keep the axial toad value within an o{}erating band and oscillating around 
a design value. This design value is ascertained from the weights of the rods, the pump, and the fluid 
column inside the tubing; while the pump minimum submersion level is taken into account. 
Minimum submersion of a pump is the lowest fluid level in a well that may be allowed for maximum 
fluid production while causing no risk to well completion- 
Diagnosing important events such as pump-off, gas-lock, over-torque, rod string fractures, 
worn pump elements, and stator swell is readily enabled by the present invention. It will be 
understood that progressive cavity pump systems are based on a production rate, which is directly 
related to the pump (rotor) speed. This speed is the same speed of the drive head shaft. 
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Accordingly, by measuring the system RPM, performing an axial load aided production 
optimization control action is possible. The system tau^t by the present invention automatically 
adjusts the RPM as required. Once the RPM adjustment is achieved, the system conducts a RPM 
check to assure that the correct adjustment was completed. The present invention contemplates that 
any of the several different ways to achieve RPM measurement including a magnetic pick-up, a serial 
connection to a variable speed drive, etc, may be used. 

By measuring the motor current, the present invention may conveniently perform several 
operating analyzes such as: predicting or detecting mechanical load conditions, stuck pumps, pumped 
solids, gas, etc Furthermore, motor current values provide the abili^ to monitor the system*s 
electrical system corresponding to balanced loads, phase loss, etc Of course, under certain 
Qioimstanoes, the system needs to be shutdown as soon as possible. Such conditions typically may 
include excessively high well head pressure, excessive torque (current), extremely low or high axial 
loads, etc Tlierefore, a means of performing an immediate, automatic shutdown is necessary. 
Alarms may be provided as pre-emptive warnings of potential system problems. 

By monitoring flow rate and pressure, it has been found to be advantageous to control the 
amount of chemical or diluent volumes introduced into the well bore or flow line. The chemicals may 
be required, as is common practice in the art, to prevent corrosion, paraffin build-ups, etc Diluents 
are used to control well fluid viscosity. These injection controls can manage the output rate of a 
chemical injection pump or the opening percentage across a chemical or diluent injection valve. To 
accomplish such controls, a means of adjusting such valve such as a valve actuator is required. The 
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instrument system conipE&s with any additional instrumentation commonly used for oil well 
applications. For exaraple, Instruments Society of America (ISA) standards are supported. 
Similarly, the present invcntioa supports other designs as well. 

The expert system (^'Master" computer system) aspect of the present invention performs real- 
time analyses heretofore unknown in the art (see FIGS. 5 and 6). Relative to axial load, down hole 
pump performance is dia^iosed. It is now feasible to measure or calculate the pump slippa^ during 
operation, and then toconpaie this slippage with either an empirical value or with a theoretically 
calculated value to dedde whether or not the pump is worn out. As will be appreciated by those 
skilled in the art, predictive and preventative maintenance may be significandy improved. The 
present invention enables the whole pump system to be protected from operating under extreme 
conditions such as over-toique, overload, and under-load. In addition, system failures such as parted 
rods, flow-line leaks, obstructions, or stator-swell are readily detected. It is also feasible to monitor 
fluid level or flow rate targets. By calculating and monitoring real- time pump intake and discharge 
pressures, and measuring fluid gradients and fluid levels, the pump speed may be automatically 
adjusted to track a parti oJar fluid level or a conrelated flow rate. It is within the teachings of the 
present invention that suchi calculations may be made either by the expert system or may be selected 
by the user. 

As will be undentood by those skilled in the art, pump rotor RPM is a cnidal control 
variable for the system contemplated by the present invention. A means of adjusting RPMs is 
necessary to control the production rate and fluid level of the oil or gas well. It is, of course, 
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important that the control method implemented complies with prevalent standard instrumentation 
communication protocols, e.g., Modbus RTU, Modbus +, TCP/IP, 4-20 mA, etc). Some of 
these devices are the variable speed drive, servomotor pulley system, etc., or servomotor pulley system 
commonly used in the art, i.e., mechanical variable speed control device wherein servomotor adjusts 
the speed output of variable pitch pulleys. Thus, the expert system disclosed changes the pump 
speed appropriately to optimize pump operating conditions. Then, the system waits for the well to 
recover, depending on the well completion dimensions, flow rates, etc During this recovery period, 
the instrument system does not provide any dynamic calculated information because the well is 
unstable and tlierefore, any calculations would be inaccurate. Nevertheless, the present invention 
continues to provide the measured information, and the basic control rules continue to operate the 
system. 

The expert system taught by the present invention provides control rules that suggest potential 
pump failure, well condition, and provides alarms and even shutdowns to prevent damage to 
expensive pump system components. It also provides explanatory messages and concomitant 
information to well operators, explaining well productivity and pump performance preferably in real- 
time. As will be hereinafter described in detail, the expert system contemplated hereunder preferably 
uses fuzzy logic to generate a unique set of operating parameters for each oilfield application, and 
collects and analyzes historical data into a knowledge base for enabling long term decisions about 
pump and well performance to be conveniently made. The present invention also exploits self- 
generated parameters to provide practitioners in the art a dynamic knowledge base having a library 
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of pump and well performnce data which is continuously updated. Thus, the expert system taught 
by the present invention affords practitioners a novel synergy wherein an already comprehensive 
analysis may be inherentlyimproved because all prior analysis and recommendations are incorporated 
therein for subsequent lilcevnse analysis and recommendations. 

As will be hepeinafcer described, this analysis and recommendation capability of the present 
invention may be interfaced with or transferred to other well applications, thereby enabling maximum 
optimization to be achieved in minimal time. Once the expert system of the present invention 
establishes rule sets fiporn a icnowledge base, then these rules may be applied to other wells. The 
expert system may be configured to interface with any current industry standard automation scheme. 
As will be appreciated by those skilled in the art, extensive data entry and analysis are provided for 
use assessing oil well equipment performance and the like. It will also be appreciated that 
embodiments of the present invention have been developed to analyze analog system values for 
several components integral to the electro-mechanical surface drive system dedicated for progressive 
cavity pumps. As will be hereinafter described in detail, mechanical relationships for these analog 
values are developed intostiitable mathematical algorithms used to generate useful values for dynamic 
pumping conditions. TTie computer system taught by the present invention processes and records 
these values, wherein the analysis is in real-time. It is an advantageous feature of the present 
invention that well adjustnnerts are made continuously. The preferred embodiment of the present 
invention presents these vadues and analyses in a user friendly human computer interface, 
programmed in Visual C-I--+ for use under the Microsoft Windows environment. Thus, it will 
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become evident to those conversant in the art» that the technology for finding these analog values, 
methods of analysis, and computer system designs are heretofore unknown to progressive cavity pump 
drive systems and the like used in oilfield production and recovery, operations- That is, practitioners 
in the art have used downhole pressure transducers connected to a surface processor or connected to 
dedicated sonic apparatus for ascertaining real-time fluid level information. However, such 
oonventional systems are not only comparatively expensive, but also afford limited capabilities because 
merely a single variable indicative of fluid level may be analyzed; sudi systems inherendy lack the 
ability to accurately diagnose dynamic pump conditions and to provide an expert system for 
optimizing pump performance as contemplated hereunder. 

Referring to FIG. 5, diere is depicted General Loop 300 of the present invention which 
fiinctions as an automatic control closed loop for the PLC. As will be hereinafter described in detail, 
the General Loop 300 comprises steps of scanning variables comprising analog inputs and 
parameters 3 1 0, then diagnose procedures 320 for optimization as contemplated hereunder, and a 
procedure for reporting variables 330 to the Master Expert System or SCADA System. Now 
refcnring to FIG. 6, there is seen the steps comprising diagnose procedure 320. More particularly, 
certain rules are observed, i.e., Rule 1 (340) and Rule 2 (350), and certain analyses are performed, 
i.e., Analysis I (360), and Analysis 2 (370), are successively executed as will be hereinafter 
described to comprehensively and dynamic assess pump operating condition and performance. 

Exploiting the logic and various functions of the expert system incorporated by the present 
invention requires that a methodology be implemented that is contrary to vAiBt has been heretofore 
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performed by operators i rtlie art. It w\\\ be appreciated that the underlying concepts and logic are 
premised on a principle ascertained from extensive laboratory and field testing: a functional 
relationship between aial load and fluid level has been found to be behaviorally accurate and 
consistent when properl^r adjusted for the cumulative effects of all the variables that impact this 
fmctional relationship. Indicative of the efficacy of the computer system taught by the present 
invention, the relabonslmip, as developed by the instrument expert system, can be and has been 
expressed as a simpler ecpialion than nnight normally be anticipated by those skilled in the art. 

Procedurally, a aeries of steps is prerequisite to properly initializing a particular oilfield 
application. During thecal Irati on of such an application, initial fluid levels are obtained from sonic 
level testing or other methods known in the art. As will be understood by those skilled in the art, at 
this time the fluid level and axial loads become known values. Then, several fluid level tests at 
different fluid levels ar^made to properly characterize the functional relationship between axial load 
and fluid level. The eaip-ett system of the present invention then generates a polynomial equation to 
represent this f^ndamenalal lelationship. It has been found that the appropriate formula may be 
effectively ascertained ihtrough conventional linear regression. Of course, it may be advantageous to 
calculate the impact on the adal load and fluid level relationship by other variables with polynomial 
equations generated dam g this calibration step. It is a feature and advantage of the present 
invention that the capad^of the instant computer system to use these self-generated polynomials and 
actual calculated value for fluid level improves the reliability of the fluid level values, per se. 

The present iii\r€nEioM) lias accomplished this daunting multifarious task of real-time data 
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collection and concomitant processing and analysis, and consequent adjustments to pump operation 
in a manner heretofore unknowrn in the art. As will be understood hiy those skilled in the art, the 
present invention exploits the ability to obtain suction pressure from the axial load supported by the 
drive head. More particularly, the total fluid level derived from axial loads supported hy the bearing 
disposed in the drive head may be represented by the following expression: 

Wr + /?p + ( P.k + Px-h ( LrGs)){Alt - Ar) - PcAr " La 

LVLs^ ^ — (1) 

ArGs 

where: 



LVL, = Fluid Level . ft. (Calculated) 

L. = Axial Load - lbs. (Measured at the load cell) 

W, = Weight of Rods in Fluid - lbs (Calculated) 

Rp = Rotor Pull down - lbs. (calculated) 

= Wellhead Pressure - psi (measured) 

P, = Friction Pressure - psi (calculated) 

L, = Length of Rods - ft. (measured) 

G, = Fluid Gradient - psi/ft. (calculated) 

Ar = Pump Rotor Area Cross section - in^ (measured) 

Ay = Rod Area Cross section - in^ (measured) 

P^ = Casing Pressure - psi (measured) 



It has been found that knowledge of an operating well's dynamic fluid level provides 

-31- 



CA 02231321 1998-03-09 



invaluable insight into pump perfomance. Indeed, dynamic fluid level enables the relationship 
between pumping rate aj^d veil productivity to be observed in real-time, wherein pump performance 
may be monitored and maacimized to engender not only optimal oil or gas production, but also well 
completion and a comprehensive knowledge of reservoir condition. It will also be appreciated that 
the fluid level value, as compaued against the pump's theoretical output capacity, additionally suggests 
the condition and performaace of the pump system. As with any pump, as a downhole pump wears 
out, operational efficiency deteriorates. In conventional dovmhole pumping operations, pump output 
is manually adjusted by appropriately varying the pump RPM in order to strive to maintain 
maximum well productivi^. It should be evident to those skilled in the art, that with no real-time 
knowledge of well conditions, analysis of well and pump performance is a time-intensive and elusive 
effort. From a practical vantage point, without having the benefit of automated real-time analysis as 
contemplated under the present invention, catastrophic reduction in well production is the only 
realistic way in which the practitioner may see changp requirements for pump production- 

GenenJIy, a well operator must use a periodic expensive and inconvenient sonic fluid level 
test to obtain tlie dynamic fluid level. However, due to the inconvenience and expense of these tests, 
the well operator usually estimates the required production rate in an attempt to provide maximum 
well productivity. Additionally, there is a serious potential of excessively low fluid levels damaging 
an oil well's productivity. This known potentially disastrous consequence understandably causes 
conservative action by an operator, thus reducing the probable maximum production of any 
particular oil well. Similajriy, as is well known in the art, pumps are used in gas wells to relieve back 
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pressure on the formation by removing fluid from the well bore. The back pressure produced by the 
fluid in the well reduces or stops the production of gas. Downhole pumps are usually placed below 
the well perforations to reduce the fluid level until the gas flows out of the well. Typically, the well 
will continue to flow until back pressure recurs due to the continuing inflow of water into the well. 

Accordingly, under the present invention, suitable design parameters must be established in 
order to properly configure the various components of the integrated computer system so that real 
time pump performance may be monitored and pump optimization adiieved under the influence of 
an automated exp>ert system. The optimization aspect of the present invention has broken down into 
two primary control units. Thus, the preferred embodiment is configured using a primary control 
system comprising the ACU. It will be appreciated that the primary function of the ACU is to 
provide a detailed evaluation of well operations, containment of the expert system, and evaluation 
of the stored data obtained from a local controller placed locally at the well site, i.e., obtained from 
a PLC or from a RTU- The PLC/RCU thus provides local control of the expert system taught 
hereunder via access to a plurality of variables that are dictated either by downloading a 
corresponding plurality of values firom the ACU or by values programmed that are manually 
programmed during system assembly and installation. 

These two control units taught by the present invention require different design parameters. 
The ACU, functioning as the "brain" of the system, interactively interfaces with the practitioner/user 
and requires significant computer programming. As will be hereinafter described, a plurality of 
objects comprise the computerized optimization system, and either perform a particular function or 
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are a prerequisite for recursive human interaction. Unlike the ACU, the PLC/RTU contemplated 
hereunder, generally requires simpler programming. Indeed, the RTU is programmed during 
assembly using and during operation using downloaded values from the ACU. 

It will appreciated that, in the preferred embodiment, the software implementing the ACU 
requires substantial user input because of the naturally required elements associated with well sites 
and in order to achieve the elaborate and robust pump monitoring and performance contemplated 
by the present invention. This plurality of input variables (and associated resulting calculations) and 
interrelationships, establish the functional basis for the underlying computer program modules and 
concomitant objects as will be hereinafter described in detail. The following enumeration of elements 
correspond to lx)th user-required inputs and calculated values. It will be understood that many of 
the user inputs are preprogrammed in the software as a user interface for data entry. 

For instance, such preprogramming of user inputs would occur if variables are known to be 
standard and acceptable in the industry. Contrariwise, other prerequisite variables mandate end- 
user input to provide design details necessary for implementing all of the capabilities contemplated 
in the automated system hereunder. Of course, it should be clear diat the following elements only 
represent the basic object definitions and means for entry into the computer programs implementing 
the preferred embodiment. In many practical applications, polynomial linear regressions are used 
in the software to establish accurate values derived from proven field tests, and are indicated with 
their elements. Utilization of these prerequisite elements is illustrated by pseudo-code hereinafter 
described. Furthermore, integration of these elements and objects, and related calculations may also 
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be readily seen in this pseudo-code. For convenience* the following elements enumerated without 
concomitant comment generally refer to user or programmer entries; elements are categorized within 
their intended objects contemplated by the preferred embodiment which is implemented in the C-h -h 
language. It will, of course, be understood that implementation may be achieved in any suitable 
programming language on any suffidendy powerful and versatile desktop or portable computers. 

The ACU design objects comprise well object, pump object, rod object, well- completion 
object, drive head object, gas anchor object, surface motor object, reservoir object, communications 
library object, tubing object, casing object, flow line object, gas separator object, load cell object, 
protocol object, surface equation object, variable frequency controller object, downhole equipment 
object, pump optimization system object, drive head object, automation surface equipment object, 
coupling object, section rods object, continuous rods object. 
Well objects elements and behavior are broken down as follows: 

well id name 

static fluid level: calculated via the static pressure: 

static bottom hole pressure: calculated via the static fluid level 

pressure of well flowing: calculated via the formula 

pwf= (perf depth*specific gravity oil*water gradient) + well head pressure 
actual producing rate: calculated via the pump operating rate; formula to be descnbed 
hereinafter 
desired rate: 
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maximum rate: r :: • crcalculated using the standard Vogel & Darcy equations, well known in the 
leum industry, for loor de evaluating inflow performance. If the productivity index, pi=0: 
I Hftftf bubble point pressure > static pressure -> Vogel 
I BEE!Zlse standing (Vogel or Darcy) 

I BEBse via the productivity index and static pressure (line equation since pi 

o Qcacurve is a straight line) 
bottom hole teirrrr=:nm!i|«peratupe 
well head tempecqB^^'i-^^^^^ 
productivity inocmdHteex: 
bubble point pncDrresssssLire 

well head pressta^ssLiULire: Measured on real time operation 

calculated via correlations from the gas separator, 
gas oil ratio: cajsjzradJt.o.culated to compensate on real time operation 
wor (water cut • * rt • cDor water oil ratio) 
oil API gravit>nc=ttty* 
gas specific gra'& — ^ai*.wvity 
water salinity 
H2S ait 

bottom hole visdciii&icrrostty 
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well head viscosity 
aromatics percentage 
sand presence 
maximum operating speed 
minimum operating speed 

dynamic fluid level : calculated via the load cell measurement on real time operation 

load c ell measurement - buoyancy factor * rods weight (air) whp + fl 
pump constant * pump area * fluid gradient fluid gradient 

(2) 

calculated with the pump intake or discharge pressures and the DeltaP. 

well load operational pressure: measured by the gauge in the load cell 

multiplied by the pumping area results the DeltaP of the pump 

friction losses on tubing: calculated by: 

A = pow((0.0292 * qO), 1 .85); (3) 
B = pow(d, 4.8655); (4) 
Tf = . 1 8 * (length of tubing or pump depth/ 1 00.0) * (a / b) ; (5) 

fluid gradient: calculated via a linear regression of practical empirical curve upon start-up 

volumetric factor: calculated via correlation 

radius of the well 
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draining horizoa tal ra^dius 

buoyancy factor dfctemnined on start-up procedure via a linear regression 
formation volumetric factor 
casing: see casin^ciefinitioa 
PCP pump object elein«nt:ssuicl behavior are broken dov^m as follows: 
brand: 
model: 

maximum rate @0[iead @ 500 rpm 
maximum head: 

maximum rate @niax Head @ 500 rpm 

rotor diameter 

eccentricity 

hp @0 head @^0Orpni 
hp @ max headi @ 500 rpm 
number of sta^ s 

operating delta jiesssiire: cadculated via the load cell measurement 

Lo ad - Buoyancy Factor * W eight of Rods 

Deltas = " — ' r — : (6) 

Constant of Pump * Pumping Area 

pump intake pr«&ure: calculated via correlations based on the calculated fluid level, i.e., 
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based upon underlying design process 

calculated via correlation based on the DeltaP, and the well head pressure 
pump setting depth: 
minimum submersion of pump 

operatiog RPMs: calculated via a parametric fionction depending on the operating DeltaP, 

and upon the desired or actual rate: 

vb=qmax5/50O; (7) 
aa= (qmax5-qapmax5)/pow(pmax,2) ; (8) 
Rrpm= ((qmax poz+aa*pow((2.308*DeltaP),2))/vb); (9) 

operating rate: calculated via delta pressure of pump and operating rpm; tested by measuring 

the presence of fluid 

elastomer: determined by conditions of the well, e.g., aromatics, sand, etc 
k constant for pumping area 
minimum starbng-up torque 
Rods object elements and behavior are broken down as follows: 
friction losses on rods: calculated via the correlation 
diameter 
roughness factor 
linear iveight densi^ 
weight on air 
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Completion ob)ect elements and behavior are broken down as follows: 
perforation depth 
total vertical deptb 
casing 
tubing 
rods 

flow lines 

true vertical depth cf every piece 
Drive head object elemenls and behavior are broken down as follows: 
brand 

model: recommended by the computer program via the maximum axial load to be generated 

by the system 

maximum axial load 

maximum hp 

minimum rpm 

maximum rpnn 

Gas anchor object elements and behavior are broken down as follows: 
efficiency 

Surface motor object elements and behavior are broken down as follows: 
brand 

-40-^ 



CA 02231321 1998-03-09 

model size: calculated by the computer program via the maximum hp to be managed by the 
drive head 
maximum hp 
maximum temperature 
maximum current 

actual temperature: measured via a temperature transmitter 

actual current under operation: measured via an amp meter in one of the phases of the motor 
balance condition: calculated via the value of the current in every phase 
Resenroir object elements and behavior are broken down as follows: 
name 
field 

geographical address 

reservoir thickness 

vertical permeability constant 

horizontal permeability constant 

number of wells 

wells 

Communications Lib. object elements and behavior are broken down as follows: 
protocol 

hardware platform (layers of the operating system model, serial port. Network port) 
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configuration 

Tubing object elements and behavior are broken down as follows: 

number of sections 

internal diameter 

outside dieuneter 

roughness factor 

weight linear density 

length of the section 
Casing object elements cuid behavior are broken down as follows: 

number of sections 

internal diameter 

roughness factor 

weight ]inear density 

length of the section 
Flow line object elements and behavior are broken down as follows: 

length 

internal diameter 

pressure drop-down: calculated by the program via the gas separator pressure and 
head pressure 

calculated having one of the above and using the correlations. 
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roughness factor 

average angle firom the horizontal ( up "4- "» down **-'*) with the vertices on the well head 
Gas separator object elements and behavior are broken down as follows: 
volume capacity 
separation pressure 
temperature 

actual pressure: measured via a pressure transmitter 
Load cell object elements and behavior are broken down as follows: 
total area to load 

pressure of hydraulic fluid: measured by the pressure transmitter on the load cell 

total axial load : calculated via the pressure of the hydraulic fluid multiplied by the total area 

to load 

maximum axial load measurable: calculated via an analysis on the load cell 
Protocol object elements and behavior are broken down as follows: 
communications configurabon 
frame 

error messages 
error cliecking 
function codes 

Surface eq. object elements and behavior are broken down as follows: 
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motor 

variable frequenoncncy controller 
Variable frequency contilrMitroller object elements and behavior are broken down as follows: 
input voltage 
output voltage 
fipequency range 9^ 
frequency step 

Downhole equipment oUo object elements and behavior are broken down as follows: 
pump 
gas anchor 

PCP Optimization Sys^''cystem object elements and behavior are broken down as follows: 
well 

automation surfalWaoe equipment 
downhole equipiq ^pment 
PCP well modeatzdel 
Drive Head object elemm^infients and behavior are broken down as follows: 
it is a drive heao&*ad but has 
load ceil 

Automation surface eqiup-quipment object elements and behavior are broken down as follows: 
variable frequen«fisncy controller 
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aim drive head 

Coupling object elements and behavior are broken down as follows: 
external diameter 
length 

friction losses on couplings: calculated via a correlation 
Section rods object elements and behavior are broken down as follows: 
friction losses on rods: calculated via the correlation 
number of sections 
diameter 
length of section 
roughness factor 
linear v/eight density 
weight on air 
couplings 

Continuous rods object elements and behavior are broken down as follows: 
friction losses on rods: calculated via the correlation 
diameter 
length of section 
roughness factor 
linear weight density 
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weight on air 

Classes relationships for a reservoir are broken down as follows: 

a reservoir has wUs 
Classes relationships fora-^U are broken down as follows: 

a well has an artificial lift equipment; a well has surface equipment and downhoie equipment 

a well has a completion; a well has casing tubing rods, couplings 
Classes relationships forartifidal lift equipment are broken down as follows: 

an A.I.M. (Drive Head with load cell) has surface equipment 

an A.I.M. has dovwnhole equipment 
Classes relationships fordovnhole equipments are broken down as follows: 

it has a pump 

it has a gas anchor 
Classes relationships forsi^ace equipment are broken down as follows: 

it has a motor 

it has a variable (re<;uiency controller 
it has a head 

Classes relationships forcompietion are broken down as follows: 
has a casing 
has a tubing 
has rods 
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has couplings 
has a flow line 

Classes relationships for PCP artificial lift equipment are broken down as follows: 

it is a kind of an a.I. equipment. 

it has a PCP pump as downhole equipment 

it has a drive head as surface equipment 
Classes relationships for pump are broken down as follows: 

has maximum rate 

has a maximum head 

has speed 

Classes relationships for PCP pump are broken down as follows: 
it is a kind of a pump 
it has a rotor 
it has stator 

Classes relationships for gas anchor are broken down as follows: 
it has an efHdency 

Classes relationships for an axial load measurement drive head are broken down as follows: 
it is a kind of a drive head 
it has a load cell 

Classes relationships for a PCP well model are broken down as follows: 
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it has basic rules 
it has expert rules 

Classes relationships for an optimization PCP system are broken down as follows: 
it has a well 

it has automation surface equipment 
it has downhole equipment 
it has a completion 
it has a PCP well model 
Classes relationships for an automation surface equipment broken down as follows: 
it is a kind of a surface equipment but 
it has a variable frequency controller 
it has an aim drive head 

Referring again to RG. 4, PLC/RTU 400 receives four analog inputs through field 
instrument manihall 10. More particularly, field instrument marshall 10 receives analog input from 
each of flowline pressure transducer 1 5, casing pressure transducer 20, motor temperature transducer 
30, and RPMs provided by variable frequency controller 60, inverter, or a RPM magnetic pickup. 
It will be understood that, when using a variable frequency controller or an inverter, the variable is 
contemplated to be provided by the Modbus serial port. It will be appreciated that all of the analog 
inputs should preferably be 4-20 mA inputs. Thus, under the present invention, the RTU will 
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simultaneously process pressure information transmitted from each of the load cell and well head, 
motor current and temperature information, and RPM information. 

The plurality of variables incorporated into the rules and formulas and algorithms taught by 
the present invention and implemented into the preferred embodiment are enumerated as follows: 
WHPH: High Rangp Value for the Well Head Pressure. 
WHPL: Low Range Value for the Well Head Pressure. 
LCPH: Hi^ Range Value for the Load Cell Pressure. 
LCPL: Low Range Value for the Load Cell Pressure. 
CALC 1 : System Axial Load 
CALC2: Operating DeltaP 
CALC3:Wei<jhtofRods 
CALC4: BOUYANCY FACTOR 
CALC5: Operating Fluid Level 
CALC6: Nominal Load 
CALC7: New RPMs. 
Kl : Weight of Rod String on Air 
K2 : Friction Losses 
K3 : FLUID GRADIENT 
K4 : PUMP Area. 

K5: The Optimum Dynamic Fluid Level. 
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IC6: Optimum Well Head Pressure 

K7: Optimum RPM. 

KB: Maximum RPMs. 

K9: Minimum RPMs. 

KIO: Pump DeluP Desired 

Kll:Qmax5 

K12:Qpmax5 

K13: Maximum Head the Pump is capable of managing. 
K14: Load Cell Area in Square inches. 
KI5: Maximum Rate at wKich the well can produce 
K16: Step for incremenhng or decrementing the RPMs 

Kl 7: Recovery Time for the well to expect the load to resume to its normal band 
KI8: Nominal Current for the Motor or VFC, etc 

As will be hereinafter described in detail, these variables are subsumed into a plurality of rules, 

analyses and p rocedures that perform the advantageous pump optimization disdoaed in the present 

invention. The convention used herein for identifying these rules, analyses, and procedures are 

enumerated as follows: 

Rule 1 : Parted Rods Detection 

Rule2: Current Detectioii 

Analysis I : Load Analysis 
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Analysis2: Well Head Pressure Analysis 
Procedure I: Pump Slippage Calculation Procedure 
ProcedureZ: Current Detection. 

ARRAY I : Array containing the Load measurement during the pump slippage procedure 

According to the preferred embodiment, new RPMs are recommended and set by the 
computerized expert system. This may also be provided via the Modbus port to a slave variable 
frequency controller CVFC"), etc As hereinbefore described, the standard of 4 to 20 mA should 
preferably be supported. The system will have four digital outputs for triggering automatic shut- 
down of the pumping assembly. It will be appreciated that normally the PLC manufacturer 
provides a minimum amount of digital outputs- It is, of course, contemplated by the present 
invention that there should be four such output in order to have sufficient room for three optional 
digital outputs that might be needed. Similarly, three digital inputs should preferably be provided 
corresponding to a "No Operation" condition, a *'Motor Overheat" alarm, and a "Current or 
Power Supply" failure alarm, h should be understood that a memory area should preferably be 
reserved for the constants to be used in this process. Most of these variables are determined during 
the start-up procedure and may be downloaded or overridden the Master upon operation. 

According to the present invention, regarding scale considerations, calculations are performed 
on the current measurement registered by the PLC (as an analog input) via the formula: 

Actual Value = C^^^ent Measuremeat(mA) ^ ^^^^ ^ 25 * IRV - 0.25 • 

1 6 
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(10) 

wherein Current Measviremenl is the analog input value; HRV is high range value; LRV is low 
range value. HRV and LRV are constants of the PLC and can be set by the Master Expert 
System of the present invention and then be downloaded to the PLC. 

System Axial or Thrust Load C'AI") is to be calculated via the formula: 

Load = LC_Pressure*LC^AREA (11) 
or Load = Strai n Gauge output 

wherein LC_Pressure is the pressure registered by a pressure transmitter (4-20 mA) used with a 
hydraulic load cell; LC^AFIEA is a constant determined by the operator of the PLC. This load 
value initially equals 3 1.42 square inches. Changes to this parameter should preferably be set by 
downloading another configuration from the Master Expert System by writing on certain Modbus 
address. Under the preferred embodiment, the system load value is reported to the Master via the 
Modbus port. 

The Pump DeltaP is determined from the formula: 

DP= (System^Load - Weight.of^Rods)/PUMP_Area ( 1 2) 

wherein DP is Pump DeltaP; System_Load is calculated as hereinbefore described; 
Weight_of_Rods = Weight_of_Rod_String^on_Air * Bouyancy_Factor; (13) 
and PUMP^Area is a constant that is downloaded by the Master Expert System. It will be 
appreciated by those slcilled in the art that Bouyancy_Factor is typically determined by the startup-up 
procedure as a function of RPM. Thus, in order to know the buoyancy factor, it is necessary to 

-52- 



CA 02231321 1998-03-09 



measure the RPM. It will be understood that there is a memory map stored in the PLC for holding 
arrays with corresponding values for both RPM and buoyancy factors. A representative memory 
map generated by the PLC of the preferred embodiment is: 



RPM 


Bucyancy Factor 


0 


0.47 


100 


0.43 


200 


0.40 


300 


0.35 



If the RPM value is in between values stored in two successive rows of this array, then a linear 
extrapolation is necessary to ascertain the buoyancy factor value. For example, if the RPM in a 
certain moment is 150, then the Buoyancy Factor is obtained via the formula: 



(0,40-0.43) 



(14) 



Then, the BF = 0.4 1 5. Under the present invention, this value is reported to the MASTER via 
the Modbus port. 

Under the present invention, fluid level is calculated from the following formula: 

(Pump DeltaP - Well Head Pressure - Friction Losses) 



Fluid Level = 



Fluid Gradient 



(15) 



wherein: Pump_DeltaP: is calculated as hereinbefore described; Well_Head_Pressure is an analog 
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input; Fluid Gradient is a constant to be set upon start-up by the Master; Friction Losses is a 
constant determined by the Master on the start-up» and may be settable during operation. 
It will be appreciated that the underlying design of the expert system taught by the present invention 
is based on some theoretical values generated from formula (1) described hereinbefore. It has been 
discovered that the present system may use a polynomial regression formula to provide accurate fluid 
level values. The boundary values for generating this polynomial regression formula is established 
during the start-up sequence in the field. During start-up, the fluid level is actually measured 
preferably using sonic equipment. The computer system taught by the present Invention evaluates 
axial load, flow line pressure, casing pressure, and known variables against sonic fluid level values 
to generate an appropriate regression formula. The result of this generated formula is compared 
against theoretical values derived from the original formula. This error checking is used to provide 
accurate fluid level values and relationships regarding other variables. Once the polynomial formula 
is thus established, accurate low-level values are provided to the Expert System and, subsequently, 
to the user in the field. 

It should be clearly understood that this model for pump design and optimization is generated 
on an individual basis for each well site. It should also be understood that this plurality of fonnulas 
and mathematical functions are programmed both in the ACU and the PLC/RTU. To provide 
accurate fluid level and concomitant information during the life span of pumping systems 
contemplated hereunder, re-calibration should be preferably performed periodically using sonic fluid 
level measuring devices. 
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Under the teachings of the present invention, nominal load constitutes a reference value for 
the system load and is used v/hen the basic rules are being executed on the PLC. 

Nominal Load = Pump Area * (Flo * Fluid Gradient + W HPo + Friction Losses) 

(16) 

+Weight of Rods 

wherein: PUMP_Area is hereinbefore described; FLo is optimum fluid level, which is a constant 
determined by the Master on start-up and settable under normal operation; Fluid_Gradient is 
hereinbefore described; WHPo is optimum well head pressure considered to be normal under 
operation, which constitutes a constant to be determined by the Master on start-up and can be 
changed during operation; Friction_Losses is hereinbefore described; Weig^t_of_Rods is 
hereinbefore described and corresponds to the buoyancy factor calculated using the value of RPMo 
to correlate; and RPMo is the optimum RPM, which is a design constant set by the Master on start 
up and settable during operation by the Master- 
Relative to the System RPM parameter, there are two constants to be set upon start-up, and 
settable during operation: MAX_RPM and MIN_RPM. As vrill be understood by those skilled 
in the art, MAX_RPM specifies the maximum RPMs allowable for the pumping system and 
MIN_RPM specifies the minimum RPMs for operation. It will be appreciated that the minimum 
RPM during normal pump operation is zero. 
The formula for the RPMs is 
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Qmax + A*(2.308 *DP)' 

RPM = 500*-^^ r z (17) 

QmaxS 



wherein: Qmax is the maximum rale possible for the pump system, or is another value settable by 
the Master^ and remains constant under normal operation; DP corresponds to 
Pump_DeltaP_Desired which is the value for the Delta? that the rules of the present invention are 
recommending to be set next — this value can be determined either by a rule or by the Master Expert 
System; QmaxS is the maximum rate the pump is capable of managing when operating at 500 RPM 
@ zero head, and is a constant determined by the Master on start-up; A is a factor determined by 
the formula: 

^ (Qmax5-Mpmax5) 

A = 2 

Pmax^ 

wherein: Qpmax5 is the maximum rate that the pump is capable of managing when operating at 500 
RPM @ maximum head, and is a constant determined by the Master on start-up; Pmax is the 
maximum head of the pump, which is a constant determined by the Master on start-up- The value 
of the RPMs under the preferred embodiment is contemplated to be set by the Modbus port as a 
new frequency value. It will be appreciated by diose skilled in the art that, for every VFC, the new 
frequency is determined by formula as a function of the requested RPM; optionally it will be an 
analog 4.20mA output. In both cases the set-point may vary depending on the %vay the Slave device 
for adjusting RPMs understands the command, whether it is a new set-point, or a certain increase 
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or decrease of the current value. If there U an increase or decrease of this current value, it is 
necessary to know the current RPMs in order to know the differential value thereof. 

The application of this plurality of objects and variahies may be conveniently illustrated using 
pseudocode. For example, in Analyrsis 1 , the load analysis aspect of the present invention considers 
every value for the load between the Hi-Normal-Limit = 0.85*Nominal_Load and The Lo- 
Normal Limit = K15*Nominal_Value to be considered normal, Henoe, for loads mthin this 
comfort 2one there is no alann generated the present invention- For events which occur outside 
of dns zone, hoivever, alarms are triggered: (NominalJValue contesponds to CALC6) 
Events: HI*Normal«Limit Exceeded: 

countl ~ counti +1 

ifcotmtl >3 Shut_Down and Reset count Generate Alarm and 

Exit Loop 

else Decrease RPMs (Decrease AOl) 

Wait for Reojveiy Time 

Lo-Nornial-Limit : 

count2 = countZ + 1 

if count2 >3 SHUT_DOWN and Reset count Generate Alarm and 

Exit Loop 

else Increase RPMs (Increase AO I) 

Wait for Recovery Time 
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Continue General PLC Loop 
It should be noted that, within this Ioop» the Increment of RPMs is a constant to be set by the Master 
upon start-up, and can be changed during normal operation. The Recovery Time contemplated 
under the present invention is a period of time during which the load will be measured, but during 
which the event loop is not going to be executed again. It should be evident that this procedure is 
observed in order to enable the pump system to recover from the new RPMs set-point. Once the 
Recovery Time has expired, of course, the event loop will be executed again if another Limit 
Violation occurred, even though the normal scan for the PLC is not interrupted by such an event. 

For normal operation, start-up of the computer system taught by the present invention is first 
configured by setting the values for every constant or parameter required by the PLC to perform the 
logic contained in the General Loop (see FIG. 5). These values include HRV, LRV, and 
Weight_of_Rod_String_on_Air which are enumerated in the following: 

HRV: High Range Value for the analog inputs. For the Load, the HRV is given by the formula: 

HRV = L5*(Rods Weight in Air * Buoyancy Factor @ ORPM + 
(Pump^ Settings Depth * Fluid. Gradient * Pump. Area)) 

LRV: Low Range Value for the analog inputs. For the Load, the LRV is given by the formula: 
LRV= (Weight_on_air_of_Drive_.Head^Shaft) (20) 
Weight_of_Rod_String_on_Air free weight of the whole rod string 

It should be noted that the HRV and the LRV for the well head pressure transmitter are settable 
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by the ^4aste^ during this start-up procedure. It is contemplated that a memory area for 
downloading all such constants will be established in the PLC. 

According to the prevent invention, a memory map should also be generated upon start-up 
to determine the buoyancy factor depending upon the RPMs. This determination is assisted by an 
Echometer or a Sonolog system in order to approximate the Fluid Level. It has been found to be 
advantageous to conduct several tests from 0 RPM to 500 RPM to arrive at a suitable Buoyancy 
Factor calculation: 

(LC - Pump Area * (FL + W ell Head^ Pressure + Friction Losses)) 

BFx = — * X RPM 

W eight. on_ Air. oL Rods 

(21) 

wherein: BFx is the Buoyancy Factor @ the RPM X; FL is the fluid level registered hy the 
Echometer or Sonolog; and the rest of the variables are already designed. Then, a memory map 
should be established in the PLC, so that the Buoyancy Factor may be determined whenever it is 
necessary for the pump performance and optimization contemplated hereunder. 

Another variable is the fluid gradient which is a constant to be determined during the well 
sizing or design phase. Based ypon an operator's practical experience, the Fluid Gradient should 
generally considered to be constant. Nevertheless, the variations on the axial load due to compressed 
free gas in the tubing fluid column should not be very high because of the low weight of the free gas. 
Experimentations in field have generally supported this theory. PUMP_Area is a variable that 
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corresponds to pumping area considering the effect of the rods. FLo is the value of the optimum 
dynamic Fluid Level which is preferably determined by knov^ng the value of the desired rate. Then, 
by using the pump performance curve shown in FIG. 12 and with the values of the desired RPMs, 
the head across the pump, i.e., the DeltaP of the pump, is calculated. Once the Pump DeltaP is 
calculated, the optimum dynamic Fluid Level may be ascertained: 

FLo= DeltaP-Well_Head_Pressure-Friction_Losses 
Fluid_Gradient 

DeltaP - Well Head^ Pressure - Friction Losses _ 

FLo = ^ : — — z (22) 

Fluid. Gradient 

Similarly. WHPo is the optimum well head pressure considered to be normal under 
operation. This pressure is typically determined by the Gas Separator and then correlating pressure 
to the well head. If the well is an old well, then WHPo might be set by the operator. 
RPMo is the optimum RPM which is determined from knowledge of the optimum FLo or the 
optimum DeltaP: 

Qmax + A*(2.308 *DP)' 

RPM = 500*-^^ ^-^—^ (23) 

Qmax5 

^/IAX_RPM is the maximum RPMs allowable for the pump system, v/hether it is limited by the 
VFC, the drive head, or by the well operator. MIN^RPM is the minimum RPMs allowable for 
the pump system. Pump_DeltaP_Desired is the DeltaP across the pump when producing the 
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desired production rate @ the desired RPM. and is determined based upon the pump curve. 
While, as will be appreciated by those skilled in the art, all of the desired RPM values are ideal, such 
values are necessary for system start-up. Variable QmaxS corresponds to the maximum rate that 
the pump is capable of managing when operating at 500 RPM @ zero head. Similarly, QpmaxS 
is the maximum rate that the pump is capable of managing \^en operating at 500 RPM @ 
maximum head. Pmax is the maximum head that the pump is capable of managing. 

It should, of course, be clearly understood that this start-up procedure should include the 
calculations of the values for each and every variable or constant described herein. These 
calculations are based upon the information provided by the well operator, maximum rate, API 
gravity, etc. A start-up related consideration is Rods Pull Down Load calculation. The Rods Pull 
Down Load is the load that is caused by the rod stretch when the pump starts running. While there 
is no known reliable way to determine the RPDL, with some practical assumptions, this value can 
be measured. Under the preferred embodiment, the Echometer Reading is considered to be the 
actual pull down load measurement. Accordingly, the calculated fluid level with the axial load will 
be corrected as appropriate. It has been generally found that the fluid level measured with the axial 
load is less than the Ediometer measurement. If, in any moment the Echometer reading is EFL, and 
the Axial Load reading is AFL, then a correction number — the RPDL — will be determined: 

RPDL= (EFL - AFL)*FLUID_.GRAD1ENT*PUMP^AREA (24) 
It has been found that, to assure that pump performance is prop>erly monitored and optimized, 
RPDL should be measured upon start-up and at a very low speed, e.g., less than 1 00 RPM, to 
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make sure that the higltrJbJiest value is taken. 

The Elxpert Swv^— vstem tierein considers fracturing of these rods via a Rods Fracture Rule. In 
particular, If System^IUE .'ILoad < Weight_of_Rod_StringLon_.Air*BF @ 0 RPM 

Shut_mZDovm and Generate Alarm and Exit Loop. 
It should be understock <==)d tbat the buoyancy factor, BF, used in diis correlation is calculated at zero 
RPM to assure that tHrl-iheie is, indeed, a fracture condition. For a Pump Slippage condition, the 
Expert System proceeQC»= :ds as follows: 
If Not Shut_DOCZZ)owi 

Shut_lIIIIZDowTi vhile System_Load>0.85*Nominal_Load or 

Time_IBt_ Elapsed < 30 seconds 

Continrm mue to Measure Load 

Store L__Zl*oad and Time Startup 

Loop 

Exit Procedure iu ■ »c 

The following i psetido-code illustrates the Expert System's procedure for current detection: 
Hi-Nora ««mnaJ-Liniit Exceeded: 

Shut_Dawi and Generate Alarm and Exit Loop 
Loop 

Lo-No>i <=DmiaJ-Limit : 

Shut^Dowi and Generate Alarm and Exit Loop 
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Loop 

Unbalanced G^ndition 

Shut_Down and Generate Alarm and Exit Loop 
Loop 

Power Supply Failure 

Shut_Down and Generate Alarm and Ejdt Loop 
Loop 

Referring now to FIG. 5, diere is shown control and optimization algorithms taught by the 
present invention. Relative to axial load, M^>* I^oad is the load value at the point of shut down due 
to overload. Min. Load is the load value at the point of shut down due to underload. DV (Design 
Value) is the load at which the well is under normal conditions and producing at its optimum rate. 
Decrease RPM corresponds to the conditions under which the system will automatically decrease 
the pump RPMs until the load value is within the normal operation band. This normally equals 
1 . 1 *DV, \^erein the customer makes the selection with recommendation by the expert system. 
Increase RPM corresponds to the conditions under which the system will automatically increase the 
pump RPMs until the load value is within the normal operation band. This normally equals 
DV/I • I , wherein the customer makes the selection with recommendation by the expert system. The 
load values are compensated due to impact of factors aflecdng hydraulic loads, such as flow-line 
pressure, casing pressure, and mechanical friction in the downhole elements. Under normal 
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operating conditions, the atxial load is to remain within its band once the system has reached the 
RPM at which the target fluid level occurs, and speed variations may take place according to load 
variations. 

Maximum Axial Load corresponds to the conditions under which the pump differential 
pressure is high due to a lo^ pump suction pressure, suggesting that the pump production rate is too 
high. RPM is decreased until the axial load goes back to normal after the recovery time. Several 
adjustments are attempteci preferably with the number of adjustments being selected by the user. If 
the system fails to recover^ a shutdown command is generated. Minimum Axial Load corresponds 
to the conditions under which the pump differential pressure is low due to a high pump suction 
pressure. Pump speed increments are issued to increase the axial load. If this fails, after several 
attempts, a shutdown is generated. If the axial load reaches a value lower than Minimum Load, 
referring to FIG. 5, theev^ent of parted rods is detected, and a shutdown conunand is immediately 
generated. Well Head Rressure, Current, Temperature Analysis corresponds to the conditions 
under which there is a shutdown command or alarm generated when any of these variables exceed 
selected maximums or mini mums. 

It has been found that there are advantageous general tips for application of the PLC Rules 
described herein. First, the Increase or Decrease RPMs command of the Expert System can be 
overridden by a Shut-Dovvn command, but the contrary is not possible. Second, the PLC will run 
the Pump Slippage Procedure upon request from the Master. Next, all of the computer system 
constants, analog inputs and outputs, digital inputs and outputs, and step-wise calculations (variables 
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calculated in the PLCs formulas or Rules and procedures) must be retrievable and settable by the 
Master Expert System or by a SCADA system, upon operator request. Furthermore, upon 
detection of any of the three digital input alarms hereinbefore described, the system taught hereunder 
will be shutdown, and the alarm will be generated for the Master to retrieve it the next time it requests 
an Information Log. 

Under the present invention, two WHP upstream of the choke and downstream analuysis 
proceeds thusly: 

Loop 

Measure all variables 

If (Well_Head_Pre8sure_Up8tream = Well_Head_Pres5ure_Downstream) 
Go ahead for the Regular Analysis 

else 

If (WelLHead_Pressure2 IS Normal) 

If (Well_Head_Pressurel is Low and Load is Hi and 
CurrentNormal) 

Free Gas across the Pump and Decrease RPM 
If (Well_Head_Pressurel is Low and Load is Hi and CurrentHi) 

Sand Across the Pump and Decrease RPM 
If (Well_Head_Pre5surel is Normal or Hi and Load is Lo and 
CurrentNormal) 
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W 



Well^Flowing Naturally, Shut_Down and Report 
If (Well_Head_Pre5surel is Normal or Hi and Load is Hi and 
CurrentNormal) 

Well is losing Fluid Level and Decrease RPMs 
herein W e 1 1 _ H e a d _ P r e s s u r e 1 corresponds to 
WELL^HEAD^PRESSURE^^DOWNSTREAM and WelLHead_Pressure2 corresponds to 
WELL.HEAD_PRESSURE.UPWNSTREAM. 

It will be understood that Recovery Time contemplated by the present invention is the time 
a well needs to respond and stabilize when it is undergoing a speed change. Thus, Recovery Time 
is defined as the time needs to displace the complete fluid column above itself to the surface flowline. 
Under the preferred embodiment, this time is calculated as follows: 
Volume of Fluid Column (Barrels) 



Recovery Time = 



Daily Rate (Barrels) 

3.14/4(Tubing ID)' * PSD 
Daily Rate (Barrels) 



, (25) 



86.400 seconds. 

Other correlations incorporated into the preferred embodiment include, based upon the 
hereinbefore enumerated variables and parameters: 

CALC1 = An*K14 (26) 
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CALa= (CALCl- CALC3)/K4 



(27) 



CALC3 = Ki*CALC4 



(28) 



For CALC4, the hereinbefore illustrated memory map will be used to determine the value for this 
variable. In actual practice, die map preferably consists of two arrays of 1 0 values each. As will be 
understood by those skilled in the art, there is a corresponding relation between each member of the 
first and the same member of the other. Depending from die value of INPUTl , diere is a different 
CALC4. That is, to know the value of CALC4, it is necessary to measure the INPUTl . If the 
value of INPUTl is in between two of the same array, then a linear interpolation is performed to 
ascertain the value for CALC4. For example. If the INPUT I at a certain moment is 1 50, then 



CALC4 is determined to be 0.43 + (0.40- 0.43)/(200-100) * (INPUTl- 100) or 0.415. 



wherein for CALC3 in this context, the Buoyancy Factor is calculated using the value of K7 instead 



CALC5 = (CALa-AI2-K2)/K3 



(29) 



CALC6 = K4*(K5*K3+K6+K2)+CALC3 



(30) 



of INPUTl. 



CALC7 = 500* 



KI5 + A M2.308*K10)^ 
Kll 



(31) 



where A corresponds to a factor determined by 



A = 



Kn>K12 
K13 



(32) 
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It will also be understood that the value of CALC7, as a new frecuency or RPM value, is to be set 
by the Modbus port or via an analog output (AO I). Under the present invention, there is a formula 
for every VFC to determine what the new frequency will be as a function of the requested RPM. 
Optionally, it will be an analog 4-20mA output. In either case, however, the set-point may vary 
depending on the way the slave device used for adjusting RPMs comprehends the command — 
whether it is a new set-point or whedier a certain increase or decrease of the current value. If the 
slave device interprets the command to change the value, it is, of course, necessary to know the 
current RPMs relative to what the differential value is being changed. 

Calculation of the friction losses in the annular space between rod string and tubing during 
operation conditions implicates several parameters including Current_Rate, HoldUp, Internal 
Tubing Diameter, Rods External Diameter, Density of Fluid, Viscosity of Fluid, Length of Section, 
Pump Setting Depth. For r=0.00 1 corresponding to the minimum value for tubing relative rugosity, 
and b corresponding to Extemal_Diameter/Intemal_Diameter, i.e., the ratio between diameters, 
the diameters correction factor for calculating the Reynolds Number (Kb) is: 



for vAndi the maximum value is 1 . The Reynolds Number Correction Factor (Z) is calculated as: 



Kb = + 



log(fe) 



(33) 



Z=(l-b)2*(l-b2)Ab 



Z = (l-b) 



Kb 
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The Hydraulic Diameter is calculated as: 

IdnternalDiameter^ - External Diameter') | 

Hd = ~ 

(Internal Diameter + External Diameter) 

then the flowing area is obtained as follows: 

^Internal Diameter^ - External Diameter^) | 

3.1 4 1 6 * ;; ^ 

4 

factor Q=Current_Rate/HoldUp, and average speed in (m/s) 

Velocity = Q 

(350.62*FlowingArca) 
wherein Effective Reynolds Number Ref is ascertained by 
406.86*Velocity*Fluid Densily*Hd/(Fluid Viscosity^z) 

if Ref > 2000, then the condition is Turbulent Flow 

S =ColeBrook_Factor (reef, rr); (Colebrook Factor) 
Accordingly, Friction Losses may be established by 

(62.4/1 2)*fr*(3281 *Velodty)2*LENGTH^OF^SECTION/(2*Hd*gravity^ 

62.4 0 Length_of_ Section 

* ff * (3.28 1 * Velocity) ' ♦ r — 

12 2 * Hd * gravity 
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else , the condition is Laminar Flow and the Friction Losses may be established by 
8 * Viscosity + q + Length, of. Section 



Internal Diameter 

51 1.8505* 10000* Kb* 3.1416 *( : ) 



(40) 



2 



Once, Friction Losses are determined, viscosity and density may be calculated. 

For the Colebrook_Factor calculation, this calculation is based upon a loop that is to be 
generated until the following condition is achieved: 

|dfj <ftol AND I y| <ytol 
where ytol = 0.0001 and ftol = 0.000001 . According to the teachings of the present invention, 
df and y are calculated every time the loop is executed. It has been found that the Loop is executed 
at most 20 times if this condition is not reached before. It will be appreciated that this procedure 
guarantees that a convergence value is reached for the Loop and the Colebrook Factor f. Loop 
Calculations are as follows: 

Relative. Rugosity + 9.28^ 

log( ' 7=-) 

Reynolds. Number * vf 1 

Y=:2* = +-?= (40) 

^ log(lO) VT-LU 
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-1 

yp = (2*f'-5) " 

9^28 

9 28 

(Reynolds. Number * f ' ^ * Relative. Rugosity + — '- — - — -777-) * log( 1 0) 

(Reynolds. Number * (f)) 



(41) 

then (if=y/yp; and the new value for f is adjusted before the loop is repeated 
f=f-df. 

If the number of times the loop has been repeated exceeds 20, the loop is exited and the Colebrook 
Factor is the last value of f. 

Relative to the calculation of viscosity, to ascertain viscosity during operational conditions, 
a loop is also generated to calculate an average value between two heists, i.e., between two points 
of the tubing string. For instance, a Well head height (0) and Pum Setting Depth may be selected 
heights. Under the present invention, the Loop divides the tubing string into smaller pieces and then 
the Average Temperature is calculated in that piece. Next, the Average Viscosity is calculated under 
these conditions and the value is added to an accumulator in the end of the Loop. It will be 
understood that when the calculations have been made across the whole string, the Accumulator is 
divided by the String Length. This loop is mathematically equivalent to calculating the integral of 
the viscosity with respect to the height: 
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-■p^dViscosityci,.,, (42) 



dh 

laitial Depth 



Thus, in this instance, • the Initial^Depth corresponds to the Well Head point and the Fmal^Depth 
corresponds to the Puirmmp Setting Depth, 
int nn; 

double h,imudh.dh,tct* t,m,niu,inuO,whdm=0,rdm=0Jnitialjremperature; 
nn=10; 

The differential of Heii:=ight is the total length divided by 1 00 { 1 00 pieces) 

dh=(Final_I>»^pth.Initial_Depth)/nn; wherein the Initial value for the differential of the' 

viscosity is 0. 

iniudh=0; 

According to the presc9«^nt invention, the initial temperature is calculated with: 

InitiaL Temperatu lore = WelL Head. Temperature + 

(B. aottom. Hole Temperature -WelL Head.Temperature) (43) 
^^P^*^ * Perforation. Depth 



muw= 1 ; Watoczer Viscosity 
Dilute Fraction in liqiuruid 

r = Dilute_Fr;TBiraction/(Dilute_Fraction+Oil_Fraclion); 
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Water Fraction In liquid 

fagua== Water_Cut/(Water_Cut+Dilute_Fraction+Oil_Fraction); 
Calculation of the viscosity of the oil with the function Viscol with a plurality of parameters: 
vc=viscol (Oil_api, Wdi^HeadJTernperature, WeH_Head_V>scDsity, Bottom_Hole_Temperature, 
Bottom_Hole_Viscosily, Initial_Teinperature). If the dilute fraction in the tubing (Bottom Hole 
Dilute Injection) is different from zero, then the viscosity of it is calculated as well: 
vd= viscol (Dilute API Gravity, Well_Head_Temperature, Well_Head_Dilute_Viscosity, 
Bottom_Hole_Temperature, Bottom_Hole_Dilute_Viscosity, Temperature). 



Auxiliary Variable = e 



r*(log(Iog(vd+{).7)))+(l.r)»(log(log(vc.0.7))) 



(44) 



And the mix viscosity is determined from: 



vm 



Auxiliary Variable 



0.7 



(45) 



The density of the mixture is : 



Dilute Density = 0.000343 * (60 - Current. Temperature) + 

141.5 



(46) 



(131.5+ Dilute API) 



141.5 



(47) 



Oil Density = 0.000343* (60 -Current. Temperature) + 



(131.5 + Oil API) 
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dm = r*Dilute DensitJis asL i: ii . .ity 4- ( l-r)*Oil_Density 



Accordingly, the Mix V 
Mix Viscosity = 
Continuing to calculatealjsat* 
muO = Mix Viscosity y^zjy 
Loop 

Calculat&lA*la 



Viscos ity is obtained: 
the Promedium Viscosity, the initial Mix Viscosity is 



.ateFempieratupe as hereabove at the depth H: 



TempC = WelLHQHDK Mnr.TTl eal. Temperature + 

Bottom. Holers A e . - Tenn|>erature -WelL Head. Temperature 



H 



lerfo ration Depth 



The Mix Viscosity is C3 ass o 
at the height H: 

mu = H = 
According to the presema 
i.e., mu and muO, 

m=(mirtTliT 
and the difference bet^»«-€tl : 
function is accumulate 5i.mxi& 



calcula-ted again at the temperature TemC which corresponds to uie vi 
Mic V^i&cosity 

sntinvention, an average value is calculated between the current two vise 
s=tvwce n the two heights being managed Initial. The argument of the i 
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imudh = (m*dh)+imudh: // integral of Mu*dh 
Then, the value of initial viscosity is adjusted to calculate the next piece of tubing string: 
muO=mu; 
End of Loop 

It will be understood that the Loop is executed from the first point, i.e., the Well Head to the Bnal 



point, i.e., the Pump Setting Depth. Next the Integral is calculated: 

Viscosity = imudh/(Final_Depth-Initial_Depth); (51) 

The procedure in the Initial Viscosity Calculation herein referenced as VISCOI. 

VISCOI: 

tkl =(Well Head Temperature - 32)/1 .8+273. 1 6 (52) 

tk2 = (Bottom Hole Temperature - 32)/l .8+273. 1 6 (53) 

dl= 0.000343*(6a.Well HeadTemperature) + l41.5/(131.5+APD (54) 

d2 = 0.000343*(60.Bottom Hole Temperature) + 141.5/(131.5 + API) (55) 

vl = Well Head Viscosity/d 1 (56) 

v2= Bottom Hole Viscosity/d2 (57) 

a = (log(Iog(v1 +0.7))4og(log(v2 +0.7)))/log(tkl/tk2) (58) 

trk=tk2 (59) 

vr=v2 (60) 

tk= (Current Temperature-32)/1 ,8 +273. 16 (61) 
Current Temperature depends from the depth 
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auxiliary variab le = c * ^^^> 

V = gt*"^'*"^^**'"''** «J - 0.7 (63) 

d = 0.000343* (fi)-Curre nt_ Temperature) + 1 4 1 .5/( 1 3 1 .5 +api) (64) 
Viscosity = ( v)*(.<l) (65) 
It will be understood thattfciis promedium viscosity is calculated assuming linear dependence between 
density and temperature. Of course, this procedure might be used instead of integrating it: the 
assumption would be that the diluent and the oil are totally misdble. If there is no dilute injection 
bottom hole, then the visioslty caui be integrated via VISCOl hereinbefore described. 

According to thepiresert invention, the following procedure has been found to be useful for^ 
calculating fluid density m dervwll operating conditions. The plurality of input parameters necessary 
are: last theoretical or pnctical punnp slippage, pressure at the point the density is to be calculated, 
temperature at the poinltlie density is to be calculated, dilute fraction at the point density is to be 
calculated (0 if not injedefld) , qqg = volume of free gas above the pump (estimated by correlation), 
oil specific gravity, dilute STpwific gravity, gas specific gravity, water specific gravity, water cut, psep. 
tsep, rssep, bottom hole static pressure, gas oil ratio, bubble point pressure, bottom hole temperature. 

Procedure: 
Calculation of Oil Volumetxic Factor 

bov= Volume tri^Fa<tor_Lib(Pressure, Temperature, Gas_Spedfic_Gravity, Oil_API, 

psep, tsep, rsscp , Bottom_Hole_Pres8ure, Gas_ Oil_Ralio, Bubble_Poinl_Pressure, 
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Bottom_Hole_Temperature) (66) 
Calculation of the Gas Volumetric Factor: 

bgas=Gas_Volumelri_Factor(Pressure, Temperature, Gas_Spedfic_Gravity) (67) 
.^ount of Liquid: 

qql = (Oil_Fraction*bov+Water_Cut+ Dilute_Fraction)/Oil_Fraction (68) 

Calculation of the Holdup: 

holdup=qql/(qqI+qqg); \^thout slippage in the ideal situation (69) 

holdup=Slippage+(l-Slippage)*holdup; with slippage (70) 

Calculation of the RS Factor for the gas: 

r8v=RS(Pres5ure, Temperature, psep, tsep, Gas_Specific_Gravity, Oil_API, rssep, 
Bottom_Hole_Pres8ure,Gas_Oil_Ratio, Bubble^Poin t_ Pressure, 

Bottom_Hole_Temperature) (71) 

Density of the oil: 

rhoo=(62.4296*OiLSpedfic_Gravity+0.076366*rsv*Gas_.Spedfic_Gravity)/(b^ 

calculated in lbs/ft3 (72) 

bl= l-Oil_Fraction+Oil__Fraction*bov; for Volumetric Factor Liquid (73) 

fw=Water_Cut/bl; for percentage of water (74) 

fd = Dilute_Fraction/bl; for percentage of dilute (75) 

fo=OiLFraction/bl; for percentage of oil (76) 
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Density of the whole liquid: 

rhol=rhoo*fo462.4296*(Water^Spedfic_Gravity*fvv+Dilute_S^^^ 

Density of the gas: 

rhog=dengasCPressure, Temperature, Gas_Spedfic_Gravity) (78) 
Based upon these predecessor parameters, the Mixture Density may be calculated: 

Density = rhol*holdup+rhog*(l -holdup) (79) 

End 

It will be appreciated by those skilled in the art, that several procedures depend from the hereinbefore 
described Density Calculation Procedure: 
Volumetric Factor 

The plurality of Volumetri^Factor Parameters consist of: Ga8_Spedfic_Gravity, Oil_API, 
psep (Separator Pressure if indicated and Well Head Pressure if not), tsep (Separator Temperature 
if Indicated and Well Head Temperature if Not), rssep (RS factor in Separator if Indicated and 
0 if not), Bottom_Hole_Pres8ure, Gas_Oil_Ratio, Bubble_Point_Pressure, Bottom_Hole_ 
Temperature). 

double sgo,co,rsv,pbt,bov,sg100; 
Calculation of the Factor for correction of the volumetric factor 
sglOO = GasSpedfic^Gravity*(l+0.1395*OiLAPI ) 

*tsep )*logl0((p8ep+14.7)/ll4.7)): 
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sg 100 = Gas. Specific^ G ravity * ( 1 +).l 595 * OiL A PI^'^'' ) * 
{tsep'-''')*loglO 



.0.466 ....... P^^P -^14.7 (80) 



114.7 



Oil Specific Gravity. 

sgo= 141 .5/( 1 3 1 .5 +OiLAPI) (81) 
Calculation of the RS Factor 

rsv=rs(Pressure, Temperature, psep, tsep, Gas_Specific_Gravify, Oil_API» rssep, 
Botlom_Hole_Pre88ure, Gas_Oil_Ralio, pb, tr) 
Calculation of the Bubble Point at the temperature: 

pbt = lTipbt(Bubble_Point_Pre85ure, Temperature, Bottom_Hole_Temperature) (82) 
if Pressure < =pbt 

// use Glasso correlation.recommended by University of Tulsa for API<20 

Gas Specific Gravity 

bov = rsv * ( — = + 0.9688 * Temperature (83) 

sgo 



bov = -6.5851 1 + 2.91329 * loglO(bov) - 
0.27683 *(logl0(bov))' 



bov = 1 + 10'"^ (85) 
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else 

b,, = * Gas_Specific.Gravity^,„, ^ ^ ^^^^ , Temperature (86) 

sgo 



bov = -6.585 1 1 + 2.9 1329 * loglO(bov) - 
0.27683* (loglO(bov))' 



CO = 



(Pressure + 14.7) 



(87) 



bov = 1 + 10^-' (88) 
co=loglO(6.8257) + .5002*iogl0(rsv) +0.36I3*loglO(Oa_APl) 

+0.7606*loglO(Temperature) (89) 

co=co-035505*logl0(5gl00) (90) 

10" 



■/ 1000/1000 (91) 



bov = bov *e (92) 

Calculation of the Factor RS: 

Plurality of parameters include: Pressure, Temperature, psep (same as hereinabove), tsep (same as 
hereinabove), Gas_Specific_Gravity, Oil_API, rsep (same as hereinabove), 
Bottom_Hole_Pressute, Gas_Oil_Ratio, Bubble Point Pressure, Bottom Hole Temperature. 
Factor for Compensating the Gas Specific Gravity: 
sglOO = Gas Specific_Oavity*(l +0. 1 395*Oil_API ) 
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*(tsep )*logl0((psep+14.7)/n4.7)): 



sglOO = Gas_Specific_ Gravity * (1+).1595 * Oil. APl"^"'") * 

psep + 14.7^ 
(tsep"^^^)*loglO( ^ ^^^^ ) 



(93) 



rss=OiLAPI/(Tetnperature-»-459.67) (94) 
Bubble Point Pressure calculation at the temperature T: 
pbt=fnpbt(Bubble Point Pressure, Temperature, Bottom Hole Temperature) 
If Bubble Point is Greater than Static Bottom Hole Pressure 
If(pbt> Bottom_Hole_Pre8Sure) 

At most Pressure can be equal to the Bottom Hole Pressure 
IfOiLAPI<=30 

rss= 0.5958 * sglOO""" *(p + 14.7)' '"»^ MO" '*"*^"" (93) 



else 

rss = 0.0315 * sglOOO.7587 * (p + 14.7)' "" * 10"^""" (96) 
If (Pressure > Bottom_Hole_Pressure) 

I I (Pressure - Bott om. Hole. Pressure) 

rsly = |Gas_ Oil. Ratio f rssep - rss |* (p^t. Bottom. Hole. Pressure) 
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Rule I for parted rods detection: 

If CALCl <K1*CALC4@INPUT1 =0 

Shut_Down and Generate Alarm and Exit Loop 
Procedure I for pump slippage calculation: 
If not Shut_Down 

Shut_Dovm 

While CALCl >0.85*CALC6 or Timc_Elapsed<30 seconds 

Continue to ARRAY[i] =CALC1 

Store ARRAY[i] and Teime StampAl [i] 

Loop 
Exit Procedure 

Thus, the system is shutdown upon request of this procedure (Procedure I). Referring to 
FIG. 7, the axial load is to be registered with its time stamp, in order to know what the load gradient 
vs time. The braking system must be functioning to keep the pump from losing the fluid column. 
During the very first 30 seconds after shutdown, the load and its time stamp is measured. As known 
the Load is given by: 

Load = PUMP_AREA*(nuid_Level*Fluid_Gradient+ WelLHead_Pressure) + 

Weight_of^Rods^on_Air*BF@ORPM) ( 1 04) 

substituting the fluid level (and the column associated to the Well Head Pressure) with an equivalent 
Fluid Column H: 
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Load=PUMP_AREA*(Hnuid_Gradieni+ Weight^of^Rods^^^^ ( 1 05) 

Now differentiating with respect the time: 

dLoad/dt = dH/dt*Fluid^Gradient*PUMP_AREA (106) 
isolating dH/dt and multiplying the annular space area O^etween Tubing and Rods: Asa): 

dH/dt*Asa = dLoad/dt* (Asa/(Fluid_Gradient*PUMP^AREA)) (107) 
the variation of the Volume V displaced due to the variation of a colunm H is dV =Asa'^dH and 
the rate associated to that volume is Q = dV/dt. Therefore: 

Qslippage =dLoad/dt*(Asa/(Fluid_Gradient*PUMP_AREA) ( 1 08) 

where dLoad/dt can be approximated to the average speed of variation (in intervals) of the measures 
taken by the PLC» upon request of this procedure. 

It will be understood that the pump slippage is to be determined under operating conditions. 
The operating condition for the pump means is when a dose value to the Nominal Load is supported 
by the pump means. On the other hand» experimentally it is considered that the well is completely 
restored 30 seconds after SHUTDOWN, These are the two main reasons why the operating pump 
slippage is measured immediately after shutdown, i.e., this is the closest value to the Nominal Load 
for the pump under shutdown conditions, but within 30 seconds after shutdown, i.e., the well is not 
restored yet. This Operating Pump Slippage may be compared to a theoretically determined 
operating Slippage (via the pump characteristics) and prematurely, a pump WORN OUT 
condition can be detected. 
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Rule 2 for current detection; 
Hi-Normal-Limit Exceeded: K18*1.30 

Shut_Down and Generate Alarm and Exit Loop 

Loop 

Lo-Normal-Limit: Kl 8*0.7 

Shut^Down and Generate Alarm and Exit Loop 
Loop 

For illustrative purposes, the following is a general step-by-step installation of a progressive 
cavity pump system contemplated by the present invention. In step 1 , the stator is attached to the 
first joint of the production tubing string. Then, the operator sequentially installs subsequent joints 
of tubing until the stator is at the required setting depth. Next, the operator secures the tubing in the 
well using conventional methods known in the art. In step 2, the rotor is attached to the first rod of 
the production sucker rod string. Then, the rotor and sudcer rod are inserted into the interior of the 
production tubing. Next, the other sucker rod sections are attached to the production rod string and 
the rotor is lowered to the stator depth. As will be appreciated by those skilled in the art. the rotor 
will pass through the interior of the stator and will then rest on the stop-pin. In step 3, die rig pulls 
the production sucker rod string upwardly until all of the slack therein is removed. Then, the 
operator marks the sucker rod string indicating its position compared with the surface tubing 
elevation. Calculations are next made to determine the amount that the sucker rods in the well will 
stretch during dynamic pump operation. The well service rig then pulls up the sucker rod string and 
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the final sucker rod element is length-adjusted with short sucker rod lengths^ i.e., with "pony rods.** 
to compensate for the expected length of sucker rod stretch and associated distance to the stator stop- 
pin. Critically spacing of the rotor in the stator has now been achieved. In step 4, the drive head 
is attadied to the production tubing and then secured to the well as required. 

It will be understood that installation of the pump optimization system taught by the present 
invention will vary somewhat for eadi individual application. The systems' remote instrumentation 
and I/O device's and concomitant processor may be housed in some various containment peripherals 
such as NEMA-rated fiberglass or metallic enclosures. Process input devices incorporated into the 
preferred embodiment as hereinbefore described consist of three transducers of various ranges. 
Additionally, the drive head contains a strain gauge. Each device will usuaDy be a two-wire 4-20Ma 
current input signal. Additional instrumentation elemjenls required based upon site. constraints or 
by the operator may be provided and connected to the system. As will be appreciated by those 
skilled in the art, local and national wiring codes should govern the installation and selection of 
shielded cables used;. Marshaling of the process input signals may be performed when a single 
multiconductor cable bridges the remote terminal unit with the process system (haying RTU, PLC, 
etc) at the well sitCi . ; ^; . ^ ... . ^ 

It will also be appreciated that the point of tap for transducers will vary from application to 
application. At the very minimum pressure, transducers should preferably be mounted in a vertical 
disposition and perpendicularly of the process line being monitored. The transducer leads, as 
hereinbefore stated, should meet local and national codes and is usually client-driven for a particular 
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application. 

To benefit (rom the teachings of the present invention, of course, the computer system 
connections including the ACU and the motor controller must be properly installed. The 
RTU/PLC should be direcdy connected to the ACU preferab^ via hard wire serial cabling or by 
using a conventional radio telemetry system; System designs for a cluster (multiple) and singular 
ajpplications may be used and, of course, vaiy significantly^ Field connections should preferably be 
made between the RTU/PLC and the motor speed control via wire analog or serial control 
conventions. 

Field testing of the present invention has enabled pump optimization to be attained to the 
extent heretofore unknown in the art. FIGS. 8- 1 1 show the results of such an actual field test. The 
values depicted represent actual responses and are not adjusted or embellished mathematically: the 
values shown represent the relationships indicated. The axial bearing loads reflect variances from 
the nominal bearing load at static well conditions. The fluid levels represent variances in feet above 
the suction of the pump compared with static fiiiid levels. Static conditions occur when the pump 
is shut down and the oil or gas well is allowed to reach a stage of stasis or equilibrium.. Test 
responses gave very good representative values for the relationship between axial bearing load and 
fluid level. 

Specifically referring to FIG. 8, there is depicted a plot of pump RPM versus measured axial 
bearing load» as RPM is increased. Thus, the relationship between measured axial bearing loads 
— corresponding to the thrust from the sucker rod and pump hydraulic loads — and pump speed. 
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As is known in the art, the pump RPM has a direct relationship to pump fluid production. As the 
RPM increases, the pump flow increases proportionately. This plot suggests an increase in an axial 
bearing load as pump RPM/flow rate increases, which proves the existence of a mechanical 
relationship to the increasing hydraulic load due to pump flow rates. 

Now referring to FIG. 9, the relationship between pump surface discharge pressure and 
pump RPM or flow rates is shown. The increase in surface pressure depicted is a result of flow line 
back pressure. It will be understood that the increased back pressure has a mechanical relationship 
to a hydraulic load generated at the pump. These results thus substantiate the requirement for 
compensation for this relationship in the prescribed calculation for establishing fluid levels. FIG. 1 0 
shows the relationship between pump lU^M/flow rate and fluid level. As the pump flow rate 
increases, the fluid level drops dearly showing the well response to increased pump discharge. 
Accordingly, this plot proves the premise that increasing or decreasing pump RPM can control fluid 
level. 

Now referring to FIG. 1 1 , there is seen the relationship between axial bearing thrust load and 
fluid level. The value for fluid level is indicated by a data line from measured fluid levels — obtained 
from sonic fluid level measurements — in the well. Also depicted therein are calculated values 
generated from the preferred embodiment as hereinbefore described in detail. Thus, the relationship 
shown proves the mechanical relationship between fluid level and axial bearing load. The calculated 
fluid level has been generated from the axial bearing load using the formula taught by the present 
invention. It should be evident that the calculated value proves the accuracy and viability of the 
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present invention for optimiung pump performance as contemplated hereunder. The ability to use 
these values from historic axial load measurements to assess* analyze, control, and predict well 
performance is heretofore unknown in the art. 

The heightened perfonmanoe achieved fay the present invention as hereinbefore described in 
detail is shown in FIG. 1 2. Performance is depicted both as a plot of production rate versus feet of 
head of water and, altemativeiy, as horsepower versus feet of head. The data shown is based upon 
100**F water. The relationship between performance and dynamic fluid level tauj^t by the present 
invention is clear. 

Other variations and modifications will, of course, become apparent from a consideration of 
the specific embodiment and illustrative examples hereinbefore described. Accordingly, it should be 
clearly understood that the present invention is not intended to be limited fay the particular disclosure, 
embodiment and examples hereinbefore described and depicted in the accompanying drawings, but 
that the concept of the present invention is to measured 1^ the scope of the appended claims herein. 
What is claimed is: 
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THE EMBODIMENTS OF THE INVENTION IN WHICH AN EXCLUSIVE 
PROPERTY OR PRIVILEGE IS CLAIMED ARE DEFINED AS FOLLOWS: 

1- In an oil or gas recovery well having pump means and 

interconnected pump control means, motor means for driving 
said pump means, drive head means, and a downhole production 
tubing assembly and a sucker rod string assembly contained 
within a well casing means, a system for monitoring real-time 
pump performance and for optimizing such pump performance, 
said system comprising: 

a plurality of sensor means disposed along said downhole 
production tubing and sucker rod string assembly for 
collecting data functionally related to said pump performance; 

a computer system interconnected with said plurality of 
sensor means and said pump control means for storing said data 
in a database, and for optimizing operation of said pump means 
by controlling said pump control means on the basis of a 
functional relationship between axial load and dynamic fluid 
level in said well . 

2. The system recited in Claim 1, wherein said 

plurality of sensor means comprises a surface discharge 
pressure sensor attached to a production line disposed 
perpendicularly of said downhole production tubing assembly 
and said sucker rod string assembly, a casing pressure sensor 
disposed on the annular space between said downhole production 
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tubing assendbly .an.d said casing means ^ an axial load measuring 
means and motor Kciurrent measuring means disposed upon said 
drive head means , 

3. The syst-em recited in Claim 1, wherein said computer 

system comprises : 

programmabl«e logic controller means electrically 
interconnected w dt-h said plurality of sensor means for 
providing local ccsntrol of said well pump conditions based 
upon a first set of rules; 

communicati oms Link means electrically interconnected 
with said programnabLe logic controller means and a variable 
frequency drive ineaxis for controlling the speed of said motor 
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means; and 

advanoed control means for providing an expert system for performing real-time 
performance and productivity analysis of said data stored in said database based upon a second 
set of rules and for providing a user interface means for enabling an operator to conveniently 
interface with said computer system. 

4. The system redted in Claim 3, wherein said computer system further comprises: 

remote control means electrically interconnected with said plurality of sensor means for 
providing remote control of said well pump conditions. 

5. The system redted in Claim 3, vAerein said second set of rules of said expert system of said 
computer system further comprises: 

procedures for pre-emptive warnings of potential pump failures and conditions for 
immediate, automatic shutdown of said well. 

6. The system redted in Claim 3, wherein said second set of rules of said expert system of said 
computer system further comprises: 

procedures for automatic adjustment of RPM of said pump means and a check-cycle for 
assuring that system RPM was properly adjusted. 

7. TTie system redted in Claim 3, wherein second set of rules of said expert system of said 
computer system further comprises: 

procedures for diagnosing well events. 

8. The system redted in Claim 3, wherein said computer system further comprises: 

procedures for generating explanatory information corresponding to real-time well 
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conditions and pump perfonnance* 

9. The system ledted in Ckiiri 3, wherein second set of rules of said expert system of said 
computer system further comprises: 

procedures for estahlishing real-time operating parameters for optimizing performance of 

said pump means. 

1 0. The system redted in Claim 3, wherein second set of rules of said expert system of said 
computer system further comprises: 

a dynamic knowiedgp base of historic well data and operational parameters compiled from 
said real-time analysis of said data stored in said database. 
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FIG. 3A 
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FIG.4 
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FIG. 5 
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FIG.8 
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FIG.9 
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FIG 10 
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HG.11 

AXIAL BEARING LOAD CHANGE RELATIVE TO FLUID LEVEL WITH CALCULATED 
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